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SYNOP«JIS 


OF 

THE THESIS ENTITLED 

"CONTPIBUTI<T*IS TO RECENT FRONTIERS IN T’lBOLOGY" 

Tribology is the name of an interdiscinlinery science which 
signifies and covers all aspects of problems related to lubrication^ 
friction and wear In its totality ^ even if parts of the problem 
belong to various disciplines of Enfineerlny and. Sciences In this 
thesis some original contributions to Tribology are embodied Ihe 
study has been mainly confined to the following aspects of lubrication 
in various trlboloyical systems 

(1) Effects of additives in fluid film lubrication 

(2) Non-Newtonian power law fluid lubrication 

(3) Kagnetohydrodynamic lubrication 

Ihe contributions have been described in ten chanters as follows 

In Chanter I, a survey of the recent contributions on the 
above three aspects has been outlined and tho prospective of the 
present work is indicated 

In Chapter II, tie generalized foynolds equation for lubricants, 
containing additives whlcli can be represented as micropolar fluids, is 
derived and its application to optimum one dimanslonal slider bearing 
has been studied by using calculus of variation It has been shown 
that the load capacity and friction force Increase and the coefficient 
of friction decreases as the parameter characterising the micro-structure 
of the base oil due to presence of additive increases 



In Chapter III tne above theory of nicropolar fluid lubrication 
has been apnlied to an externally pressurised bearinsr by using variational 
approach Hie optimun film thichnoss profile is obtained for maximum load 
daoacity of the bearing Tlie of foots of additives on the various 
characteristics of this ootinun bearing lave been investigated 

In Chapter IV the effects of viscosity variation, due to change 
in concentration of tne additive In the base oil on the characteristics 
of slider and externally nressurised bearinps have been studied It has 
been shovm that both tie load capacity end friction force increase as 
the concentration of the additive In the base oil incre&ses 

Cliapter V deals with the characteristics of non*Newtonian newer 
la\i lubricants in an externally pressurised conical step bearing and in a 
hydrostatic sten seal It has been pointed out that the step in the ^ilm 
thicloiess contributes to the increase in the load capacity of the above 
systems The effect of flow behaviour index of the power law fluid on 
the characteristics of those tribological systems is also studied 

In Otapter VI the thermal effects in externally pressurised 
bearings using a non-Newtonian power law lubricant have been studied by 
considering the consistency of the lubricants as a function of temperature 
and pressure It has been snoim that the load canacity of the bearings 
decreases as the perametors characterising the dtange in consistency increase 
Further, in the case of externally pressurised circular bearing the load 
capacity increases as the flow behaviour index decreases, but the load 
capacity in the case of rectilinear thrust plates does not depend upon 
this index The thermal effects on the flow flux of the bearing are also 
investigated 



In Chapter VII, the characteristics of non-Newtonian power law 
lubricants in the case of externally pressurised porous thrust bearing 
have been treated It is found that the load capacity increases as the 
flow behaviour index or the parametoi characterising the porosity of 
the bearing increases 

Chapter VI I J is devoted to the stud/ of squeeze film between 
rectangular and circular plates using power law lubricants where the effects 
of inertia, pressure and step in film thickness have been considered It 
has been shown that the effect of inertia and pressure is to increase the 
load capacity of the bearing, while the effect of step in the film thickness 
is to decrease the load capacity Similar effects have been found on the 
squeezing time in both the above cases 

In the last two chapters certain tribological aspects of magneto- 
hydrodynamic lubrication are studied in externally pressurised bearings 

In Chapter IX, the characteristics of an ('xtexnally pressurisod 
bearing using conducting inconpressible lubricant in the presence of 
non-uniform axial magnetic field have been investigated by using variational 
approach It has been shown that for maximum load canacity of the bearing 
both the film thickness and the magnetic field profiles should be step 
functions Hiis maximum load Increases as the strength of the applied 
magnetic field increases 

In Chapter X, the magnetohydrodynaraic extemally pressurised 
porous thrust bearings have been studied in the presence of uniformly 
applied transverse magnetic field It has been shown that the load 
capacity increases as the siqjply pressure, viscous permeability of the 
strength of the applied magnetic field increases The effects of squeezing 
have also been studied and it has been pointed out that the load capacity 
increases as the squeeze velocity increases and this Increase is enhanced 
as the strength of the applied magnetic field increases 



LIST OF SYMBOLS 


— in Chapter IX otherwise it refers to the width of the 
^1 slider bearing 

strength of the magnetic field 

volume concentration of the suspension 

volume concentration of the base oil 

dimensionless coefficient of friction 

specific heat at constant volume 

diffusion coefficient 

e}q)onential function 

Weierstrass excess function 

friction force at the slider for maximum load capacity 
film thickness 
minimum film thickness 

film thickness in regions I and II respectively 

maximum film thickness corresponding to optimum profile 

hj/hj in Chapter V otherwise thickness of the porous surface 

Modified Bessel functions of order zero of first and second 
kind respectively 

modified Bessel function of order one of first kind 

rate of mass transfer at the surfaces (in Chapter IV) 

thermal conductivity (in Chapter VI) 

step location 

rate of chemical reaction 

length of the slider bearing 

consistency index of the power law fluid 

constant consistency index 



torque in hydrostatic bearing (in Chapter IV) 

Hartmann number « (in Chapters IX and X) 

flow behaviour index of the power law fluid 

pressure in the thin film of the lubricant 

exit or ambient pressure 

externally applied pressure 

inlet pressure 

maximum uressure 

pressure at the step 

pressures in region I and II respectively 
pressure in the porous matrix 
Flow rate or flux 

volume flux \fith uniform film thickness, or constant consistency 
index 

volume flux in the porous capillary 

flux in regicm I and II respectively 

volume fluxe in x and z directions respectively 

radial variable 

outer radius 

inlet radius 

step position (in Chapter III) 

’^e 

— in Chapter III otherwise it refers to outer radius 
^i 

inlet radius (in Chapter IV *) 
average radius of the porous capillary 
squeezing time 

frictional force at the lower plate in slider bearing (in 
Chapter IV) 
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T 


temperature Ciii Chapter VI) 
u velocity in x or r direction 

U velocity of tho slider 

V velocity in y direction 

v’ velocit) in the porous capillary 

V squeeze velocity 

W, W, W! U U W_ dimensionless load capacities as defined in the 

analysis 

X = ~ in Chapters III and IX otherwise it refers to a variable 

^i in Cartesian coordinate system 

Xjj « ^ (in Chapter III) 

y co-ordinate normal to the plane of the film 

V B ~ dimensionless film thickness 

h^ 



z rectangular coordinate 

0^, Y material constants (in Chapters II and III) 
a semi vertical angle of the cone (in Chanter V) 

Y * P Cy (in Chapter VI) 

5 (qierator for variation 

w viscosity 

Pq viscosity of the base oil 

y. viscosity of the liquid dispersed in the lubricant forming 

^ emulsion (in Chapter IV), material constant (in Chapter; II end III) 

Vi 

— an index to the micros true ture character of the micropolar 
^ model 

P fluid density 

a ccaiductlvity 

TwiT* shear stress in the lubricant and porous capillary respectively 
xy' xy 

n angular velocity of the rotor (in Chapter IV) 


(0 


nacrorotation vector 



CHAPTER - I 


lETRODUCTION 


1 1 whkt is TRIBOIOGI ? 

In general, most of the lubricated systems can be considered to 
consist of two relatively moviag siirfaces (plane or curved, loaded or 
unloaded) with a thin film of external material (lubricant) between them 
The presence of such a thin film between these surfaces not only helps 
them to support oonslderable load but also minimises friction 

Since the days of Osborne Reynolds [1686 ] , various mathematical 
modelB have been proposed to study the oharaoterletics of such syetmns 
under various simplifications, Oops [1949] , Bowson [1962] and 
Cameron [1966] However, to cope with the complexities of modern 
man-aud-maohlne systems which arise in space travel, nuclear reactors 
etc , the need for new frontiers of research Involvlag lubrication, 
friction and wear has arisen To study such problems in a more realistic 
and comprehensive manner, the knowledge of the fallowing le desired 

(1) Hature of the surfaces, such as roughness, elasticity, Ihexmal 
couductlvily, hardness, pore si 1y, and its affinity such as adsorption, 
absorption, chemical reaction, etc with the lubricant 

( 2 ) Nature of the lubricant solid, liquid or gas, Newtonian or non- 
Newtonian behaviour, 'variation of vlsooslty or consistency with 
temperature, preesuro or opno entratlon and Its affinity with the surface 



(?) Effects of euxroundlog eavlxomenti external and Internal forces 
etc 

It IS obvaous from the above that to have the total picture of 
the problem} the jEnowledge of verxous disciplines of Science and 
Engineering is required In particular} If In the lubricated system} 
adsorption of the lubricant is taking place on the metallic surface due 
to chemical reaction} then the knowledge of Chemical Ihglneering} 

Chemistry, Metallurgy and Ihyelcs etc , is involved Perth sr, if the 

surfaces are rough, then to study their topography due to randomness of 
roughness, the knowledge of Statistics, Production Eoglneerxng and Metal 
Biysios etc , is required If the system forme a part of htman body- 
such as a human joint, the knowledge of Biology, Medlolne and Ehgineering 
is needed* !Do make a comprehensive mathmnatioal model of the system, 
the help of a mathematician is also desirable Ihus to have a realistic 
Investigation of euoh problems, a combined effort by solen-tlsts and 
engineers specialized in various aspects of the problem is needed Such 
an attempt did not come forth -till the middle of thle century, and an 
effort was needed -tx> generate such Interactions in the field of lubrication 
friction and wear Keeping -this in -view, the Minis -tzy of ^ohnolo^. 
United Kingdom, appointed a committee known as Jost Committee, in the 
late sixtees, to investigate the short oomlngs In the area of research 
In the above fields and Its eventual application JUx the industry It 
IS this committee vAilch, in fact, coined the new name, ISIB0IO6Y, 
whlob signifies and covers all aspects of the problem related to 
lubrication, friction and wear In Its totall-l^, even If parts of the 
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pmblem belong to many usual diaoxplineS| suoh as Ueohaolcal Sngiaseriog} 
Industrial JBhgineerlngy Chsmlstxyy Gbenlcal Eaglneering, fbyelcs; 
Metallurgy) Textile) Leather) Medical sciences etc Shus) tribology Is 
meant to be truly interdisclplln axy In nature and Is mainly problem 
oriented in the real sense of the term ) Shvtkla [1972] 

In this thesis an attempt iA made to study the oharaoterlstios 
of certain tribological systems by using thin film lubrication theory 
iOie following three aspects have been Investigated 

(1) Effects of additives In fluid film lubrication 

( 2 ) Non-Newtonian Power law fliiid lubrication 

( 3 ) Ifegnetohydrodynamlo lubrication 

1 2 EKffBCTS OP ADDITIVES IN PEDID HIM HIBEIIOaTION 

It Is needless to emphasize the importance of additives In 
lubricants beoatuse of the diverse and highly stringent environment 
(suoh as lU) mlsslleS) air crafts and other space vehicles) to which 
they are subjected Consequently) a speolallzed knowledge le required 
for careful selection of additives for each Individual application) 
fixazler [1965 ] In partioular) to improve the characteristics of 
the base oil) the additives can be used as rust Inhlbltams (Amine 
ldios|diates)) corrosion inhibitors (suljhurlsed olefines)) fire resistant 
(halogenated hydrocarbons)) detergents (Calolun/Barlum Sulphonate)) 
vlsooelty Improvers (polymethacrylate) powders of graphite and moly- 
bdenum disulfide)) E I additive (sulphurised fats)) Mblyneux [1967 ] 

Even additives in gases have also been used and the process Is known as 
aerosol lubrication) Carr et al [1965] 
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Several experimental investigations have been carried out by 
using different additives in the base oil and the reduction in wear 
rate and coefficient of friction have been pointed out, IfficCabe [1965 ] > 
Stock [1966] , Mblyneux [1967] t Bowe [1970] , Baartz [1971] In 
particular, iCalivaldle Spalvins [1971] has carried out an experiment 
and has shown that thin sputtered molybdenum disulfide film on highly 
poliedied metal surface gives low average coefficient of friction and 
long wear life at a particular load Gansheimer [1972] has reoently 
conducted an experiment on tie lubricating properties of the mixture 
of mineral oil with certain Inorganic compounds and has pointed out that 
the mineral oil suspensions of iron and zinc pyrojAoephates have hl^ 
load carrying ability and good autlwear properties 

However, very little attention has been given to develops a 
mathematical model predicting the effects of additives in l^drodynamio 
lubrication where the viscosity of the suspension has been taken into 
account Herzig [1959] has given a mathematical relation between the 
coeffloient of friction end film thickness In the case of foil bearing 
when molybdenum disulfide has been added to the lubricant Hhis 
relation has been obtained by comparing the theoretical results in the 
case of foil bearing with the experimental data obtained by adding 
additive to the oil It has been shown that the coeffloient of friction 
depends upon a factor which is a function of ooncmitraticn and size of 
the particles 

Since, viscosity is the most in^ortant factor governing hydro- 
dynamic lubrication, the effects of additives in lubrication may. In 
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general, be etudied theoretically by a nathematloal model involving 
viscosity of l^e base oil, concentration of the additive, nature of 
the additive (solid or liquid), geomstxy of the particles of the 
addxtive in the case of solids, their interaction with the base oxl etc 
Slnstein [1906] has been the first who has given a malheinatlcal foxmula 
for the apparent vxsooslty of Ihe fluid containing additive in terms 
of its concentratxon for dilute suspension He has pointed out that 
as the 0 oncentration of the additive increases, Ihe viscosity of the 
suspension increases Since then a number of generalizations of the 
above relatloa has been proposed in the case of suepensxons and emulsions 
for dilute or low concentrations, Rutger [1962] , Sherman [1962] 

Also, Seshadri and Sutera [1968] have pointed out that the concentration 
of the suspended particles change when suspension flows and the 
concentration Is governed by the equation of mass transfer and depends 
upon convection, ohemioal reaction, size of the particles, geometry of 
the system etc Ihus, in most general form, the apparent vlsooslty 
of the suspension can be represented as a series involving powers of 
concentration of the additive which may change with ohemioal reaction, 
size of the particles and the geometry of the system eto 

Further, In the oase of physical interaction of the additives 
with the base oil, such as mloro'^tatlon of the solid partlolea, the 
flow and lubrication behaviour of the suspension may be studied by 
considering them as mloiopolar fluids. Bring en [1969] , Hudimoto and 
Tolcuoka [ 1969 ] , Allen and Kline [1971] i Datta [1972] , Agarwal et al 
11972] , Khader and Vachon [1973] In addition, in the case of 
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polymer addltivesy if the oonoeutration increasea from a certain value^ 
the suapenslon no longer bebavea as a Newtonian fluid but followa the 
non-Newtonian character [see next section] 

In view of thiS) we have considered, in Chapters II and III of 
this thesis, the effects of solid particles in the lubricant by 
characterizing it as a mioropolar fluid In chapter II, the generalised 
Reynolds equation for mioropolar lubricants has been derived and the 
optimum one dimensional slider bearing la studied It is found that 
both the maximum load capacity and the corresponding frictional force 
increase as the parameter, characterizing the micro -structure of the 
base oil due to the presence of additive, increases, but the coefficient 
of friction decreases Similarly, in Chapter III, the case of an 
externally pressurised bearing with mioropolar lubricant has been studied 
by using the techniques of calculus of variation and it has been pointed 
that the maximum load oapaclly increases as the step height ratio or the 
parameter characterizing the mlorostructvire of the suspension increases 
Ihe effects of these parameters on the flow flux are also studied 

In Chapter IV, the effects of additives have been considered by 
assuming the apparent viscosity of the suspension as a linear function 
of concentration of additives by taking mass transfer parameter into 
account Ihe cases of slider and externally pressurised bearings have 
been studied and it has been shorn th*at the load capacity and frictional 
force increase as the concen-cratlon of the additive in the base oil or 
the mass transfer parameter inor eases 
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1 3 UGN-UEVraoUlMJ SOW® lAW HWH) LUBRICAOIOH 

In the case of notion of a fluid) vihen the shear stress Is 
proportional to the rate of shear) the fluid Is said to follow 
Newtonian hypothesis and is termed as Newtonian fluid lEhose fluids 

which do not follow this Newtonian behaviour are called non-Newtonian 
Strictly speaking most of the fluids express non-Newtonian behaviour 
when subjected to precise measurement) Plnkus Mathematically) 

the non-Newtonian fluids could be characterized by the following 
constitutive relation: 

Shear stress s consistency x rate of shear 

where consisteocy caU) in gaieral) be a function of both shear stress 
and shear rate Sependlng upon this relation various types of non- 
Newtonian models such as BLngham plastic ) vlsooelastiO) power law eto , 
have been proposed) Bird et al [i960] ) longwell [1966 ] Non-Newtonian 
behaviour is almost invariably observed in solutions or melts of hl^- 
molecular-^ei^t polymeric materials and suspensions of solids in llq^uid 
Some of the examples of non-Newtonian fluids are: suspension of chalk) 
paintS) printing ink) greaseS) methacrylate in water, detergent slurries 
etc ) Bird et al [i960] , longwell [1966] 

In the last two decades or so, various investigators have studied 
the lubrication oharacteristios of these fluids in various bearings , 
finkus [1960] In particular) the characteristics of Bingham plastic 
lubricants (greases) have been studied in the case of slider bearing, 

Milne [1953, 1957] > squeeze films and hydrostatic thrust bearing, 
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Osterls et al [19^5] and Batra [19731 > roller beazangs^ Saaakl et al 
I 1962 ] Kauzlarich and Greenwood [197^ have considered a generalized 
Blngbam model known as Herschel BiUklt^ Model and have studied the 
concept of elastohydrodynamlc lubrication for greases 2he MHD^thxust 
bearing using Blngbam model has also been studied by Shukla and Prasad [1967] 

Hhe viscoelastic fluids as lubricants have been studied by Milne 
fl957a3 f Selby [1958 ], Horowitz and Steidler [1960-61 ] t Wright and 
Crouse [1965 ] j Oiosh and Mishra [1968 ] , !Eanner [1970 ] f Hung and 
Muster [l970 ])ln various bearings Chow and Salbel [19711 have studied 
analytically the film thickness and pressure distribution for a heavily 
loaded line contact of rollers using viscoelastic lubricants and have 
pointed out a reduction in the flhn thickness from its Mewtonian value 

Ihe behaviour of Sower law lubidcants has been studied by several 
investigators Hg and Salbel [196^ have investigated the use of a 
iDodified model of the pseudoplastic power law lubricant in the case of 
an inclined slider bearing and have shown that the load capacity Is less 
than that of Hewtonian lubricapit Hsu and Salbel [1965] have studied 
the case of a slider bearing by using a cubic model and have pointed out 
that the load capacity and friction forc«^ are reduced but the flow rate is 
comparatively increased in comparison to Newtonian case Tanner [1963» 

1964» 1965 ] has studied the non-Newtonian lubrication theory and has 
applied this to a short journal bearing by using power law lubrioante 
The use of Bower law fluids In squeeze films and externally pressurised 
conical bearings have also beai studied, Shukla [1963a, 1964’~b, 1964“^ ] 
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In the case of conical step bearing, It bas been pointed out that the 
load capacity increases as the flow behaviour Index decreases Shukla 
and Brakasb [1969 1 have also studied the use of power law fluid in the 
case of rheostatic thrust bearing and have shown that, fbr the ataxlmvon 
load capacity of ihe bearing, the flhn thichness should be a step 
function 

The characteristics of nonoNewtonlan power law lubricants In 
journal bearings have been Investigated theoretically by Tanner [1964, 
1965 3 and Hsu [l966] and experimentally by Dubois et al [ I960] 

Hsu [1966 ] has obtained solution applicable to both pseudoplastic and 
dilatant non-Newtonian fluids, which are valid over a wide range of 
shear stress 

In view of the above, it may be noted that the effects of step in 
film thickness, consistascy variatloa with temperature and pressure, 
porosity, elasticity, roughness and inertia, in various bearings have 
not been consideri^d for non-Newtonian lubricants, Dowson et al [1966] , 
Castelli et al [1966] , Christensen [ 1969~7ll etc In this thesis, some 
of the above uffects have been studied in exttsmaily pressurised and 
squeese film bearings for power law lubricants only 

In Chaptdr V, the eff«sots of step in film thickness on the veirlous 
characteristics of externally pressurised conical bearing and hyd3?ostatlo 
step seal have been investigated It is shown that load capaoity of 
these bearings Increase with the increase in etep height It Is also 
seen that in the case of a hydrostatlo step seal, the load oapaoily 
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increases as the flow behavioui index of the fluid increases 

In Chapter VI, the effect of temperature on the consistency 
index of the lubricant la considered in the case of hydrostatic thrt:iet 
hearings (rectilinear and circular) By considering an exponential 
variation of coneistbacy with taspdratura and pressure, it is shown 
that the load capacity decreases as the flow behaviour index or the 
parameter characterizing thb consist cjacy variation increases 

In Chapter VII, the characteristics of externally pressurised 
porous thrust bearing with power law lubricants havu been studied and 
it has be.en shown that the load capacity increases as the flow behaviour 
index or the parameter ofaaracterlzing the porosity Inoreases 

In Chapter VIII, thb charaot eristic e of squeeze film bearings 
with power law lubrLoants are studied by taking inertia and the linear 
variation of oonslstency with pressure into account It has been shown 
that the effect of inertia and pressure is to increase the load oapacity 
of the bearing when the film thickness is constant lUrtber, in the absence 
of inertia and pressure effect, it has been shorn that the load oapaolty 
and time of approach decrease as the step height increases 

1 4 MAGNBTOHVDROIfmAMIC nJBRICAITlOISr 

Mognetohydrodynamic or hydromagnetio lubrioation is the study 
of lubrication process of elbctnaally conducting lubricants in the 
presence of electric or/and magnetic fields Because of high temperature 
applications such as in nuclear reactors and space vehicles, the use of 
liquid metals as lubrloants has been suggested Since the liquid 
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metals are lees viscous than oil aad are good electrical oonductorS) 
its lubricating property can be Improved by applying external 
electro-magnetic fields in the system In fact) vihen a conducting 
fluid flows in a bearing and a magnetic field is applied in a transverse 
direction to motion) a body force is generated which modifies the 
lubrication characteristics of the system In rec nt yearS) several 
scientists and engineers have contributed in the field of hydromagnetic 
lubrication by Jceeping their objective around tbe possibility of using 
different magnetic and electric field geometries to get an increase in 
the pressurization of liquid metal bearings 

In the early sixteea, Eucks and Uhlenbush [1962] and Snyder [1962] 
have studied the cboracteristios of MEO} slider bearings and have indicated 
that the load capacity oan be Inoreaaed by using conducting lubricant 
in the presence of externally applied magnetic field Since then 
several investigators have investigated the characteristlos of MHD slider 
bearings (finite or infinite) in the presence of transverse and tagential 
magnetic fields , Elco and Hughes [1962], Hughes [1963] > Shukla [1964a] , 
Ohawla 1 19661 , Prak sh [1967] ,IiaiDanalah [1966] , Kapur [1969] and 
Agarwal [1970] A generalized Reynolds equation for hydromagnotie 
lubrication has also been derived by Shukla [1963] In addition, the 
calculus of variation approach has been used to study the optimum MHD 
slider bearing by Osterle and Young [1962 ] , Kusma [1965 J , Shijkla [1970] , 
and it has been pointed out that for maximum load capaoi-^ of the bearing 
the film thickness should be a step function Hhe characteristics of 
MHD slider bearing with non unifOxm magnetio field have been discussed 
by Bodfciewicz and Anwar [1 972 > 1972a] 
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Several theoretijcal and experimental investigations have been 
conducted to study the characteristics of MHD journal bearings in the 
presence of radial as well as axial magnetic fields^ Hughes and Blco 
[1962a] f Eusma [1963f 1964b] » Shvarts [1966] j Kamiyama [1969a] j 
Kamiyama and Sato [ 1971 ] Budsinsky et al [ 1968 ] have presented an 
analysis and experimental results fOr a MHB partial journal bearing and 
have ^own that there exists an optimum Hartmann nvuaber at which the load 
capacity attains a maximum for a given applied current 

The problem of MHD squeeze film has b^en studied by Kuzma st al 
[ 1964 ] , Kuzma [1964a] , Shukla [1965] , Shukla and Erased [1965, 1967] , 
Shukla and Frukash [1966] , and the increase in load capacity and time 
of approach due to increase in loagnetic field have been pointed out 

The oharaoteristics of externally pressurised MHD bearings have 
been studied by Ihi^es and Elco [1 962 ], Shukla and frasad [1966] , 

Kriegar et al [1967 ] , Maki et al [1966,1967 ] , Kapur [1968] , Chow [1969] 
Shukla [1965a] has studied the MHD externally pressurised bearing with 
variable film tbickness and has shown that the load capacity is greater 
than that of a bearing having a constant film thickness The effects 
of inertia In MHD hydrostatic thxust bearing have also bei.n investigated , 
Bamanalah [1967] , Shukla and Kapur [1967] Kamiyama [1969, 197®, 1971] , 
Kamiyama and Sato [1972a, 1972b] , in a series of papers, have studied 
the inertia effects and influence of wall conductance on the characteristic a 


of MHD hydrostatic thrust bearings 
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In view of tho above, in Chapter IX, the char ict is tics of an 
extcrn^ly pressurised bearing using conducting incompressible lubricant 
in the presence of a non-uniform axial magn-wtic field have been investi- 
gated by using variational api-roach It has b en shown that, for 
maximum load capacity of the bearing, both th film thickness and the 
magnetic field profiles should be step functions It has also been 
indicated tnat this mcximun load increases as the strength of the magnetic 
field, applied in the region of minimuia film thickness, increases 

In Chapter X, the magnetohydrodynamic externally pressurised 
porous thrust bearings have bten studied in the presence of uniformly 
applied transverse magnetic field It has been shown that the load 
capacity inoreases as the supply pressure, viscous permeability or the 
strength of the applied magnetic field Increases !£be effects of 
squeezing have also been studied and it has been pointed d out that the 
load capacity increases as the squeeze velocity increases and this 
increase is enhanced as the strength of the applied magnetic field increases 



OHAPTEH - II 


GENERALISia) RETNOIiDS EQUATION lOR MICROPOLAR lUBRICANTS 
AND ITS application TO OPTIMUM ONE DIMENSIONAL SLIDER 
BEARING EPPECTS OP ADDITIVES IN lUBRIGATION 

The use of additives, in the form of solid particles, in the 
base oil to Improve its lubrication oharacteristlcs In bearings, is 
well known However, very little attention has been given to 
preuict these oharacteristlcs by considering the mathematioal model 
which takes into account the presence of additives In the base oil 
Einstein [1906] has been the first who pointed out that viscosity, 
the most important characteristic of the fluid, Increases as the 
concentration of the additive increases Por low concentration of 
the additive the suspension b^baves as a Newtonian fluid, but for 
very high concentration of the additive the suspension no lon^r 
follows Newtonian hypothesis and It can only be characterleed by a 
non-Newtonian model In addition, the solid particles present in the 
fluid, apart from moving with the usual velocity of the lubricant 
might have mioxo*^otation In the fluid elements Such suspensions are 
represented by lacropolar fluids Eringen [1965] , Allen and Kline 
[1971] , Khader and Vachon [1975] In this chapter, this particular 
model of suspension Is being considered to study the effects of 
additives in lubrication The generalised Reynolds equation is 
derived and the oharacteristlcs of optimum one dimensional slider 
bearing are studied by using techniques of oaloulus of variation , 
Rayleigh [1916] , Maday [1968] and Shukla [1970] 
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2 1 gmesalissd refolds equation 

The basic eqviatloas governing the flow of a fluid containing 
additives represented as mlcropolar fluid are given by Erlngen [1966] 

- V ^ CV?xV) + ~7CV^)] 

= Cl+2vi+yj^)7V V (y+Pj) 7 X V X V + V X (i> 7 d (2 1 ) 

P3 0 ) = (o^ + + y) “ yV ^ Sf ^ Zi + ^ V 2yj^ 01 (2 2) 

7 V = 0, (2 5) 

vdiere body forces and body moments are neglecteu The various symbols 
have been defined in the nomenclature and axe not described here 

Let us consider now the flow of the lubricant In a thin film 
enclosed by two relatively moving surfaces The film thickness Is 
negligibly small in comparison to the other dimensions of the surfaces} 
so that the usual hydrodynamic lubrication assumptions are applicable , 
Plhkus and Sternlloht [1960] lEie physical situation is described 
In figure no (2 1 ) 

To derive the Reynolds equation for this lubrioant^let us 
oonsldeX} as usual} the following components of velocity and 


micro -rotation vectors 




rii> ; 1 new a MirnopoiAR fliiip in a thin film 



Fin?? ONE PIMFNSIONAL SLIDFR BEARING WITH MICROPOLAH 
LUBRK ANTS 



17 


V = (uj U2 U3) 

u = (Wj 0 0)3) 


(2 4 ) 


where the compoueats of velocity and micro -rotation vectors are 
assumed small In the dlreotlon normal to the bearing surfaces 
[Ug << u^] Using the vusual assumptions of lubrication theory 

[neglecting spin inertia etc ] » from equations (2 1) to (2 4) we 
get the following equations governing the flow of the lubricEuat in 
the film 


Ci»Vi) 


(ji+yi) 


am. 


8 y 

a^w. 


a^uj^ 

+ P- 

80.3 

, 

3\2 


ay 

ax 





3 y 2 


ay 

az 


8u, 




1 

ay 


to« = 

3 


0 > 


au, 




21.J 1.J . f 


(2 5 ) 


(2 6 ) 


(2 7) 


(2 8 ) 


(Ebe boundary conditions for the above system of equations are 
chosen as follows) 


"5 

“1 “ 
w. = 


= u 

at 

y = o» 

u^ = 0 at y - h 

(2 9 ) 

= 0 

at 

y =* 0, 

y “ h 

(2 10) 

= 0 

at 

y =1 0 , 

y =h 

(2 11) 

= 0 

at 

y = 0 , 

y =h 

(2 12) 
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Solving eciuatione(2 5 ) and (2 7) and using the boundary conditions 
(2 9) and (2 12) we obtain the following expressions for and 


Uj = U + Aj[3 sinh Xy + 2y] + 'I, [$ (cosl Ay 1)] + — , (2 I3) 


Wj « ^j^Ccosh Ay - 1) + sinh Ay - 


y -HL 

2 u.yj 


(2 H) 


where = 


[ 


h 

2(2y+yi) 


^ U sinh Ah 

*1“ ‘ ‘ ■ i t mm 

^ {pCcosh Ah 1) + h Sinh Ah} 


B 


1 


1 

inh Ah 


[-JL- 

2y+yi 


IE 

3X 


A^Ccosh Ah - 1}] 1 


2 (2y+yj)yj 

Cy+yi)Y 



2 

A 


Sirailarly from equations (2 6) and (2 8) with boundary conditions 
(2 10) and (2 11) we get the expressions for and to^ as follows 


Ug = A2tS sinh A> + ?y] + Ej 3 (cosh \y 


1) + 


■XL,. jE 

2y+yi 3z 


(2 15) 


Wj^ ® A2 (cosh Ay 


1} + B2 inh Ay 


y ^ 

2v+yi 


(2 16) 


where 
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JE 


2i2v+Vi) 32 ’ 


= 


^ ^ - A2Ccosh Xh 1)] 


sinhXh 32 


Ihe flow flux per 

unit 

width in the x and 

z directions are defined 

by the following equations 


Sc “ 

h 

/ 

0 

u^dy 

(2 17) 


h 

I 

0 

u^dy 

(2 18) 


which on using eguations (2 13) and (2 15) sivo 

^3 


0„ = 


Uh 


iE cc^) 

a* ' '•i. •' > 


X 2 6(2n+yj) 3X Ml, 


(2 19) 


Q ^ Q ^ G(— ), 

® OC 2 W + Vi ) V ’ 


(2 20 ) 


where CCh) = C(' 


1, ,3 3v, h 1 a.fi cotn a. h 

r—i - h + A ^ X 

u * orn_i.li ^ o > 


2 Cv+yj) 


“ --; — = and H. = -77- 
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On integrating the equation of continuity (2 3) we have 

It ^ “1 ^ "3 ^ 

and using equations (2 17) to (2 20) we get the generalised form of 
Reynolds equation for micxopolar lubricant as follows 


dx 




— coo S + — [ — 

6(?y+vii) ax-* 82 6 (?u+mi) 


GOO 


- II5L 
32^ ■ '> ax 


(2 22) 


It is this equation vblch determines the pressure in the fUm for a 
given function h 

1 h coth a, h r? 

Since t i — L 

as 0, oj '► 0, the function 6(h) tends to h^ Then for 

= 0, equation (2 22) reduces to the following forai 


/liL- ^ {111-. s Id H 

ax ^i?p ax"^ 87 ^i2v az^ ax 


(2 23) 


This IS the wall known Reynolds equation for ITewtonian fluids 

2 2 ORTljmUM ONE DIMENSIONAL SLIDEB BEARIN6 

In this section, we study the oharaoterlstics of micro polar 
lubricants In the oase of one dimensional slider bearing whose 
configuration is shown in figure no (2 2), Allen and Kline [1971) y 
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Balram and Shastiy [ 1 972 ] | Datta [ 1 972 ] 


In this case, the generalised Beynolds equation (2 22) can be 
written ae 


dx 


< 2 — G(h) 5E-) = i 

6(2VI+V;l) dx 2 dx 


(2 24-) 


The flow flux q( = 0^) is from equation (2 19) as 



6(2y+uj) dx 


^rCh) . 


(2 25) 


which le a constant as can be seen from the equation of continuity 

Brom equatlons(2 24) and (2 25) we have the equation detemlning 
the pressure as follows: 


d£ ^ h 5 

dx ~ G(h) GKh) ’ 


where 


5 =B- . S = 


Oh ’ 




3UL(2y+1i-| ) 


f X 


X 

L 


and 


h - ^ 

h 


Integrating equation (2 26 ) and using tht. conditions) 


(2 26 ) 


p = 0 at X s= 0 

p = 0 at X = 1 


(2 27) 


we get the expressions for pressure and flow flux as fbllows 
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/ J^dx 

0 CCh) GOO 1 1 

/ dx 

0 TOO 


f' / 


r(h) 


dx = / 


_h 

-(h) 


dx 


(2 28) 


(2 29) 


Ghe load cai-acity of the bearing (per unit viddth) is given »y 

L 

W a I pdX 
0 

which can be written as 


V/ » 6 / p dx « 6 / X % dx (2 30) 

0 0 dx 

where Vf « • • ' g -’■■■' ■' 

UL‘^(2y+y^ ) 

In the following, the load capacity W is nexxmized with respect 
to the film thickness 

h >1 (2 31) 

by using the approach of calculus of variation 

2 3 ilBTHOD-I VaRIATIONAI APEROAOH WITH HIM THICKNESS AS BODNJM) OONfROI. 
VARIABIE 


Ttollowing Maday [1968] and Shukla [1970] equations (2 26) and 


(2 31 ) are expressed as 
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1 


dx G(h) G(h) 


(2 32) 



=0 , 


(2 33) 


where i/) is a real function of x When iji = 0^ h = 1 otherwise 
for non-zero t{> , h > 1 

It la required to £Lnd out the variabl s p and h such that 
W given by equation (2 30 ) attains a maxiinum 


Let us define the augmented function F as follows 


r = fp + XjCx) 4»j(x) + X2(x) <fr2(x) , (2 34) 

where X's are undetezmined Lagrange's multipliers » which are 
continuous In tht. Interval 0 £1» except possibly at a finite 

number of points 

Using equations (2 32) and (2 33 )» the function P can be re- 
written as follows 

r = 6p + X,(x) Ar'+ A-} + XoCx) {\ 1 /} (2 35) 

^ dx G(h) rch) 

2 - 

Slnoe and do not oocw in the governing equations (2 32) to 

(2 35), ^ and h are needed to be only piecewise continuous in the 
dx 

interval 0 <^x Bbllowing Leltmann [1962] , the Euler-IftgrangeB 

equations can be written from equation (2 35 ) as follows 
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d A, 

6 = , (2 36 ) 

dx 

-i {Q f + r - h r'3 + A, * 0 , (2 37) 

G 

^2 ' I ' “ ^ 38 ) 


wh<»re d&sh denotes differentiation with respect to h As the eugnented 
function I* does not Involve x explicitely,we get the following 
fundamental relation from equation (2 33) 

6o X- ^ » C , (2 39 ) 

^ dx 

where C is an integration constants 

iSie Veierstr&ss-flrclmann oonditionej which apply at a corner of 

the extremal arc wh«-re and h are discontinuous, give 

dx 


X 


1 


x*x„ 



X*X +9 

o 


(2 40) 


(6n - X, — } 
dx 


X»x„ 


* (6p - X, ^1 
0 dx f5b»Xo+0 


(2 41 ) 


where x = x^ is a point of discontinuity for the extremal and x^- 0, 

0 denote conditions Imaediately before and after x^ respectively 
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A aecessary condition fbr the existence of local maximum for the 
load capacity W la given by Weierstrass E-function, LeitUKuan [1962] , 


which is equivalent to Hie req\,ixement that 


E = 


(h-Q) 

G 


(2 42) 


must be mlaimum with respect to h 

Kcom equations (2 40) and (2 41 ) wc. see that X, and 6p - X, 

i ^dx 

are continuous at the corner x = x where and h are 

° dx 

discontinuous Since p is continuous at x = x^, from equations 
(2 40) and (2 41 ) we must have =» 0 Then integrating equation 

(2 56) we get 


A^(x) = 6(x - x^) (2 43) 

It is clear from equation (2 43) that 5. 0 in 0 <.x iLx^ 
and — ^ — Also from equation (2 37) we have 

^2 = 0 at X = X|j Now from equation (2 38) W6 I1SVY& &jLi;h6x* ^2^ ^ 
i|i 5s 0 or both are zero Erom equation (2 33) it can be seen that 
the condition ij> « 0 infixes that h - 1 for both zero and non-zero 
Xg and the condition \|> ?^ 0 gives h > 1 and Xg ” ^ for Xg* 0 

and 0 , we have from equation (2 57) 




g(V) 


(2 44) 


where ^ is the optimum value of h In this region On substituting 
this value of Q in equation (2 29) it can be noted that 1 Hence 
Xgca 0 and 0 imply that ’ll 0 Again if Xg / 0, must be aero 
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In order that the equation (2 38) may be satlsfieli and this corresponds 
to the region where h = 1 ahis suggests that h is step function 
and it changes its value where changes sign (ie at x = x^) 

It now suf fines to take the following pmflle for h corresponding to 
maximum of W 


5 = iL , 0 tx ±x^ , 

m ’ 0 ^ 

5 = 1 » ^ 5, "* » 

where h^ is given by equation (2 44) 


(2 45) 

The qualitative behaviour of 


^ 1 » ^2 ^ aeeri from figure no (2 3) 


It can also be noted from equation (2 42) that for maximum of W, 

the condition ^ = 0 is equivalent to equation (2 44) 
dh 

Using the piofilfe given by (2 45) in equation (2 29) the flow flux 
of the lubricant is given by 




Q * 


w' 


(2 46) 


^0^(1) +0 - G(^) 

Similarly the load capacity is given by using equations (2 30)| (2 44)> 
(2 45) and (2 46) as 


3(1 ^^);„o-x„)(vi) 
i^GCl) + (1 - x^) G(h^) 


<2 47) 


where W = (l + 
m 




Here x^ is still 


a ) »- *1, = “■"“I “IT 

puL ^ _ 

undetermined, eind if we maximize W with respect to x^ as well, than 

9W 

from the oondltlon — ~ = 0, we obtain 

8^ 



h 



Xq X 

(c) 


r!( ?Hn)OfTlMUM PROFll F FOR FILM THiCKNFSS RATIO 

IFAflVf Ot-OAVPUn Of A, 

( )M / I tlAflV Bl HAViOUR Of- 
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2 G(h ) 

J _ 

1-^0 " g(1 ) 


(2 48 ) 


This equation gives the step location fbr maximum load Further fron 
equations (2 44) and (2 46)fOn eliminatiog Q, we have 


- g<\) 

1^0 


{2 49) 


The equations (2 48) and (2 49) determine and h^ for maxuauai 
hllinina ting x^ from these equations we obtain the final 
equation determining h^ as 


G(1) G(h^) - {(h^-1) G'(hj^) ~ e(h^)}^ =: 0 (2 50) 


The equation (2 50 ) is solved numerically for h fbr different 
and and then the corresponding x^ ie obtained from equation (2 4S) 

The pair of x^ and h^^ when substituted in equation (2 47 ) gives maximum 
The force of friction F at the slider (per unit width) is given by 


- F 


XqL 


1.1 V • 


m 




dx 




which on using equations (2 Ij) [with u] , (2 26) and (2 45) gives 


“I'm" 


F Vi 

mo „ . 1 


Cl + T^i) {Xp Fj + (1 x„) F,} 


0-' 


y UL 


(2 51) 
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where r, - 




2a. 


JSL 


^ OO 

2a, fi tanl (a i ) 

^ " l+yj 1 “ 


and = 


' t[S„ 1) 


2a. 


2 CCD " ni 


'•’Or,) 


2a, 


1+Pi 


tanh a 


1 


AlsOftho dunensionleas coefficient of friction is given by 



(2 52) 


2 4 MhTHQD-II M/dXIMIZATION OP THE lUNCTIONAI W 

In this section the functional W Is to be loaxiuiised with 
respect ■to the film thickness ratio h > 1 by following the techniques 
of calculus of variation, Bayleigh [1918] 

Trom equations (2 26) or (2 28) and (2 30 ) the functional W 
can be rc-wrltten as, 


i; r SJL dx / dx] (2 53) 

0 GCh) 0 C(U) 

where Q la given by the equation (2 29) 

Por an admissible 'vuxiation 6h in h there would be a 
corresponding change 6W in W, the expression for vdiich can be 
written from equation (2 55) as 
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1 “ 

X dx 


= / 


0 c^ch) 


<^(h) {x 


/ 

0 r(h) 

0 G(h) 


-) 6h dx , 


(2 54) 


whwe 


»|)(h) = hG’(5) - G(h) - QOUii) 


(2 55) 


Itor maximum V the profile for h(x) should be such as to make 6W 
always negative or zero It is ae n from equation (2 54) that the 

1 xdx 
0 ^) 

sign of 6W is given by the product of 4 (h), x - ‘"U"' and 6h 

/ 

0 &(fi) 

To make 6W negative or zero, it suffioee to choose the following 
profile for h(x) 


such that 




h *= 1 f X ^ * i ■' » 


HS„) = V'(\) - 8(\) - Q S'(h^) - 0 


(2 56 ) 

(2 55)’ 


at the upper step 0 ^x ix^ and x^ is so chosen as to make the 
following xjtprtsaion positive 1 « 


* / 1 ^, 0 fori> X 

^ 0 r(h) 0 GOO “ 


(2 57) 


The beha-\lour of tCh) , given by equation (2 55)»with respect to h 
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la shown in figure no (2 4) It is noted from this graph that «f(h) 
vanishes for h = and it is negative when h<i^,(x^;^x£l) and 
positive for h > ( h^ > 1 ) Ihus ror the region 0 ^ x x^, 

^ ^ ® Therefore, SW vanishes whatever may be the 

product of the other two tenas But in the region x^^x;^l,h = 1, 
4(h) is negative and 5h is positive and therefore <SV/ is negative 
because of the inequality (2 57) 

On using the optiiBum profile in the inequalily (2 57) we get the 
oondition detenaining x^ as followst 


\ _ °<^) 

~a(i) 


(2 58) 


Similarly from the equation (2 53) 'the step height ratio 
for optimum W is given by 




°<v) 

e'(£J 


(2 59) 


Since the conditions (2 58)ani(2 59)ar the same as equations (2 48) 
and (2 1-4) respectively, obtained in the previous Section 2 5 , 
therefore this method of optLinlzation of the functional V would give 
the same results fbr various characteristics of the slider bearing as 
already obtained in the previous section 


2 5 HBSOIffS AJJB DIS3USSI0NS 

aihe variations of Q and are shown respectively in 
figures (2 5) and (2 6) It is seen that Q decreases and 
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increases as increases It is further suen that Q Increases 
and decreases as Increases for various values of 
corresponding -values of flux and load capacity fbr the classical 
Newtonian case are obtained by letting y^ 0 or -*■ «> 

Ghe variation of F Is manifested through figure no (2 7) 

HI 

Here also It Is seen that Increases as Increases The 

variation of the dimensionless coefficient of friction c^ Is shom In 
figure no (2 8) It la not*sd that c^ decreases as Increases 

and thus the Increase In f Ictlon force Is more -than compunBated by 
Increase In load capacity 




Hi 




MO ; / VPHlAnON Of frn M PIFFERENT Mj 




FIG ? 8 VARIATION OF c, WITH p, FOR DIFFERENT M, 


CHATTER ~ III 


EXTBBITAlir THBSSURISEJ) OPKMDM BEASINC WITH IWICHOTOIAH 
PUJID AS A IDBRICAHT 

In Chapter we nave studied th« characteristics of one 
dimensional optimum slider bearing by oharacterizing the lubricant 
suspension as a micropolar fluid In this chapter also) we study 
the behaviour of the lubricant auspenaioii containing additives in 
the case of an externally pressurised bearing by considering the same 
model The approach of oaloulus of variation Is applied to study 
the optimum film thiclmess profile for the maximum load capacity of 
the bearing) Walker and Osterle [1961] 

3 1 BASIC EQUATIONS 

Consider the flow of a lubrlceint suspension) characterised as 
micropolar fluid) in an externally pressurised circular bearing 
Because of externally applied supply pressure) the euepension flows 
only in the radial direction and there is no component of the velocity 
in the S^dlrectlon due to symmetry The micro-rotation vector Is 
assumed to have only e‘^(»aponent)the z-component being negligibly 
small because of small film thickness in comparison to the other 
dimensions of ths bearing The geometrical configuration Is Illustrated 
in figure no (3 1 ) ®ius by considering 


(3 1) 


V = (u, 0, 0) 

M = (o, til) o) 




FIG 3 1 EXTERNALLY PRESSURISED BEARING WITH MICROPOLAR 
LUBRICANT 
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and using the usual lubrloatlon assumptions in equations (2 l)) (2 2) 
and (2 3) we get the following equations governing the flow of the 
lubricantt 


2 

f > 3 u dp 

(P + Vi) V, — - 

1 3 j ,2 1 32 dr 


2 

d 0) 3u n 


The boundazy conditions for u and 


a 0 


0 


tt at the surfaces are 


(3 2 ) 

(5 3) 


u a 0 at z a 0 and z = h f (3 4) 

w = 0 at a = 0 and z = h (3 5) 

Sol-ving equatlon8(3 2) and (3 3) and using the conditions (3 4) and 

(3 5) we get the expressions fen? u and u as follows 


u = A[P sinh X 2 + 2z] + B(5(cosh Xz 


1 ) + 


dp 

2 V+V 2 dr ’ 


(3 6 ) 


u a ACcosh Xz ~ 1) + B sinh Xz - 


z dp 
2V+Vi * 


(3 7 ) 


where 


h dp 
2(2p+yj) dr ’ 

h dp 1 * cosh Xh 

2(2V+Vi) dr sinh Xh 



2 C2v+Vi) Vi 
X a . ■ 


Cv+VpY 
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IQie flow flux: is defined by 

h 

Q = / 2‘ffrudy 

0 

whJnh on using equation (3 6) gives 


IT h 


Q=- 


“ C(h), 


3(2uti.j) <>r 


where 


t «3 3v, h 1 

C(h) « fc' ) = u + i 

% 2Cl+yp 

r h - 2 

V '*1 “ “T" 


ttj^h coth 


or 


M. 


i VPjJlP 

1 + yj 


il - & Y = Y/P^rt 

1 4 y 0 


(3 8) 


(3 9) 


It can be seen from the equation of continuity (23) that Q 
is a constant 

Integrating equation (3 9) end using the boundazy conditions 
p a at r ts r^ 


(3 10) 


P = Pa “ ""e 

we get the final expressions for Q and p as 


P ^ 

0 / r* * (Pi “■ Pa^ 

1 xfCh) 3(2y+yj) 


(3 11) 
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X dx 
{ xG(h) 



1 xC(h) 


(3 12) 


where X = — 

iChe load oapacll^ of the bearing le defined as 


7 /^O 

V = tt(Pj^ *• Pjj) rj + / 27 rrCp p^) dr 

r 

1 


which on using equation (3 12) gives 


h 13) 


W » 


W 


R 

I “TTr dx 
1 G(b) 

R 


’'^^1“ ^a^ / dx 

1 xG(S) 


(3 U) 


where R = — > 1 and p^- is prescribed 


3 2 RRTJRtlXRaTION OR OPimUM FIIU TmCOTESS ROR MAXIMM lOiO} 
OAPAOIir 


Rrom equation (3 14 )« It can be noted that the load oapaolty 
Is a function of h(x) In this section we detennine the optimum 
profile for the film thickness h corresponding to maximum load 
capacity of the bearing 

Rrom equation (3 14) the load capacity can be re-written as 
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(3 15) 

where 

\ G(h) 

(3 16) 

and 



If 

(3 16) 

NoWy any change 6h in h would produce a corresponding change 

6W in W 

Jkom equations (3 15) the expression for 6W can be 

written as 

- ffie - pdf 

61J = 

f2 

(3 17) 

which on using equation (3 16 ) gives 



^ f F(x) 6h dx , 

r 1 G^Ch) 

(3 18) 

where 

r(x) » xf 

(3 19) 


Hure we are latereeted ia detemlnljig that profile for h(x) which 

^ iMi 

would, make negative for all possible vaxiatloas in h 

It oan be noted from the ligures (3 2) and (3 3) that <}(£) 
and G-'(h) ere positive and Increase as h inoreases fbr h ^ 1 
Otoerefore from equation (3 1B) it is clear that the sign of 6W depends 
upon the product of f(x) and 6h It can be seen that B'(x) is 



•7 
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positive for values of x close to unity and negative for values 
approaching B , Rayleigh [ 1918 ] Suppose the value of x at 
which £'(x) chaag«.s sign Is x^^ tiien from equation (3 19) we have 

*0=1 (3 20) 


ajie condition that 5W is negative fljr any adansslble variation 6h 
would be satisfied if we choose the following profile fbr film 
-chic knees ratio | 


h = h_ , 1 < X < x^ 

m ' — — o 

h = 1 , X < X < B 

’ o~ — 


(3 21 ) 


where £ Is the ratio of the maximum to ihe minimum allowable film 
thickness 


It Is easily noted that for 1 <. fih can only be negative 

and since i*(x) is positive in this range, fiW is negative Again, 
for x^j< x^ R, dh can only be positive and since f(x) is 
positive in this range, 6W is negative Ihus, it is shown that dW 
is negative througjiout the region 1 possible variations 

of h !Qie load capacity fbr this type of film thickness profile oan 
now be written by using equations (3 15) end (3 20 ) as 

\ (3 22 ) 

The relation between and h^^^ oan be obtained by using equations 
(3 20) and (3 21 ) as 
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V Pi 

t o JC X ^ f X dx 

,2 1 Glh„) x„ 0(1) 

0 

1^0 dx ^ dx 

1 xG(hn) x^ xr(l) 


which on Integration and euapUflOcitlon gives 


(3 23) 


, 

'■CD ° (B^ xj) - 2xJ to S_ ■ I ’ 

0 O V 

where and g(i) are values of G(h) when h * 

respeotively 


(3 24) 


and h = 1 


1!he equata.on (3 24) gives the relation between the step 
location and the step height ratio^for maximum load oapaolty 

of the bearing The variations of / &(1) versus pj are 

shown in figure no (3 2) It Is seen that the function S^= G(hjj^)/6(1) 
increases as or h^ increases ibr fixed !Di8 variations 

X 

of with ^ are shown in figure no (3 4)^ from which it may be 

noted that the function S. inoreases as x^ increases for fixed B 

1 o 

Since S^ increases as or ^ increases, it may be inferred 

that Increasee as y^ or ^ increases Bow from equation (3 22) 

2 

it is noted that the maximum load capacity is propertional to and 
oonsequently this maximum load increases as or Inoreases 
It is further noticed from figure no (3 2) that <J(^) / ®(*^) lecreases 
as increases^ and so with the reasoning given as above it may be 




Fib A VARIATION OF 5, WITH 


FOR DiFFERtNT R 
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concluded that the maxinum load capasity decreases as the parameter 

increases It may also he note! that the behairxour of W tends to 

m 

that of Newtonian luhricoit when p^-»■Oor^^'»■“■ It is farther 
remarked that since incrcaa«-B in«?efinitoly as a preaches R, 
the optimum load capacity increases as the step losition mervos to the 
outer red in 8 of the bearing 


From equation (j 11 ) the expression for the flow flux can be 


rewritten as 


2 ^ 1 xff(h) 


(3 25) 


Making use of equations (3 21 ) in equation (3 25 )» the expression for 
Q can be written as follows 



(1 +-^) 



(3 26 ) 


As pointed out earlier^ 1 ^ and ('(h^) increase as h^ incrt^ases and 
G(h^) > G( 1 ) for ^ 1 [see figure ni (3 2 )) IChen from equation 
(3 26 ) it may be concluded that Q increases as inorupscs for fixed 



j^sin, since 6 ( 1 ) and G(h^)) h^ > 1 y decrease ea increases 
or decreases [see figures (3 5 ) (3 6 )] , the integral on 

the right hand side of equation (3 25 ) increases as increases 
Or decreases Hence it may be pointed out that Q decreases 
as K] increases or decroases 
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3 3 GOSrCIHSIONS 

In this chapter, the characteristic s of an externally pressurised 
optimum bearing with micropolar lubricant have been investigated foll- 
owing conclusions may be drawn from the above analysis 

(1) lihe maximum load capacily Increases as the step height ratio 

h^i or the paxaiiiieter or both increase Ibis maximum load 

capacity decreases as the parameter increases The Newtonian 
case is derived by letting o or « 

(2) lor prescribed e^plled pressure, the flow flux increases as 
b^ increases and it deorsases as increases or as decreases 



CHAITER - IV 


BFjjIBTS op VISCOSITlf VaRlAlIOiT ON TH« CHARACa?JRIS!I?ICS 
OP SIIDiH AND hydrostatic BEARINGS DDE OX) OHiiNGB IN 
CONOENOEATION OP ADDITIVES 


In the previous two chapters} we have studied the oharacteristlcs 
of slider aijd extemally pressurised bearings with lubrlceinte containing 
additives by considering thein as mloropolar fluids In this chapter} 
we consider the effects of additives In lubricant through the cheucge of 
viscosity of the base oil due to presence of additive In It 

It la well knovm that the vLsooslty of the lubricant changes 
with concentration of additives whloh In turn changes with the pro- 
cesses such as diffusion} convection} chemical reaction eto present 
In the system. } Bralthwalte [1967] As early as 1906} Jineteln 
derived the following formula governing the viscosity of a dilute 
suspension of rigid spherical particles uniformly dispersed In a 
Newtonian llquxd 

V « (1 + 2 5 C^) (4 1) 

Here is the volume of the diapers ea material per unit volume of 
the base oil Since then various models have been proposed governing 
the vlsoosltles of dlSpersioue} Rutgers [1962] and emulslons} 

Sherman [1962] In general} most of these models can be written as 

(4 2) 


U - u^il ♦ XCo ♦ XjC^ * XjC’ ♦ ] 
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where A’e are parameters depending upon size, shape and nature of 
particles of the additives In particular » at low ooncentration, 
when the particles are assumed to be deformable sjdieres equation (4 2) 
reduces to the following form as obtained by Taylor [see ButgerSy 
1962] 

li = y ri + \c ] 
o 

? 

* tK 

where A * ? 5 - ■ (4 3) 

^1 

It has also been pointed out that the concentration of the 
suspended particles changes when suspension flows Seshadrl and 
Sutera [1968 ] Keeping this In view one may generalise equation 
(4 3) as follows 

U = * \C] (4 4) 

where C Is the volume concentration governed by the equation of 
mass transfer and depends upon convection y oheioLoal reaction etc 
present in the system This equation may also be assumed to be true 
when impurities are mixed in the lubricant during the operation of 
machines 

In this chapterya generalised theory of lubrloatlon Is 
formvilated which includes the effects of addltlve8/idipuriT;les 
present in the base lubricant during actual working conditions 
In partloulary the effects of addltlves/lmpurltlea on the 
characteristics of slider and hydrostatic bearings a 3 pe studied 



4 1 Basic equations 


In genei^al, most of the lubricated systems can be considered 
to consist of two relatively moving surfaces with a tnln film of 
lubricant between them Q^e preosure generated In the film for a 
lubricant with variable viscosity Is given by the generalised form 
of Reynolds equation ^ which la the outcome of equations of motion 
and the equation of continuity) tinder thin film approximation 
theory) Bowson [19621 This equation can be written as 



♦ I- tp 

ax az 


ap 

az 


0 « u 


dx 


(4 5) 


vhere 


^ J 

n=/ 21. dy r/Z.,^y 

r ^ 0 ^ 



(4 6) 


and V is given by equation (4 4) The volume conoentratlon of 
the additive Is governed by the equation of mass transfer and Is 
written as follows by using thin film approxl^oation theory, Bird 
et al [1960] 


where k' Is the rate of loss or gain in 0 due to ohemloal 
reaction) absorption) adhesion, corrosion etc She Veloolttes u 
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and w are given as follows Dowson [1962] 



Xhus to study the effects of additives and impurities in 
lubricant, the system of equations (4 4-) to (4 7) should be solved 
under suitable boundary conditions for p and C As equations 
( 45 ) and ( 47 ) are Inter-related and are highly non-linear In nature, 
only approximate solutions are possible In the following we solve 
these equations in shnplifled form to study the obaraoteristlos of 
slider and externally pressurised hydrostatic bearing 

4 2 SLIJIS B£ABI]!ra 

Let us consider the case of an infinite slider bearing as 
illustrated in figure No (4 1 ) It is assumed that Ihere is no 
ohemioal reaction in the film but there is a possibility of greater 
mass transfer at the plate y = h by processes such as ohemioal 
reaction, adsorption or adhesion At the plate y = 0, the 
concentration of the additlve/iu^url'ty Is assumed to be prescribed 
93iey may be either due to additl've present in the lubricant or due 
to wearing out of a certain solid lubricant coating on the surface 
9to simplify the matter further we assume that the diffusion due to 




f'16 t, 1 SLIDER BEARING WITH MASS TRANSFER AT y=h 
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coavectlon and conceatratloa changes with x are email In such 

a case, the eguatlon (47) governing the conceatratloa reduces to 
the following equation 



with the following boundary conditions 


(4 9) 


0=0 at y = 0 
0 

(hi - h ) 

_ do . „ . . , 1 o 

D^ = kO at y = b = h^ ^ X 


(4 10 ) 


Solving equation (4 9) and using conditions (4 10) we get a linear 
concentration profile as follows 


If there is no mass transfer at y = h, 1 e when k » 0, then the 
system behaves as if the bearing is lubricated with an additive of 
uaifonm volume concentration G 

o 

after C is determined we Can write the modified viscosity 
of the base lubricant by using equaticn (4 4) as follows 


u . n + xo^ (' )i u 12) 

Vilh this expression for u , the Beyuolds equation (4 5) in this 


case is written as 



59 


dx F 



(4 13) 


h 

where the flow flux: Q = ^ u dy = a coaataat 

Solving equation (4 13) and using the usual boundary oondltion 
p a 0 at X s 0 and p s 0 at x » L| we get the expressions for pressure 
and flow flux as follows t 



Q “ 



(4 15 ) 


Also the expressions for load oapaoity and friotional fozoe at the 
moving plate are given respeotlvely as follows by using equations 


(4 14) and (4 8) 




T «s 



(4 16) 


(4 17) 


Sino e the oono entration of additives or laipurities are fairly small f 
the functions V and G in above integrals may be approximated as follows 
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P = 


12li^ 


[1 





kh 

D-kh 


)3 


G = 


h 

2 


[1 



kh ^ 
D-kh-* 


With these expressions for P and G, the expressions for flow 
rate) load oapacity and frictional force are ^Iven as follows 


a+1 


XCo Oj. 


r s 


6 In a 

Ca-ID^ 


12 

a^l 


* V 

o C 




^ an a 


- — ; — * T- 

a+1 a-1 0 C 


where 


a =^,K 


h. kh^ „ Wh 

i.K = o « -T-^. w = — ^ 


2Jh 


h_ »“ D > “ Dh^’ ^ P UJ 


n s tt r oKA ctK ^ _ od, 

C a+1 1 ^3 a+i'^ o+l'^ 


« « {(a-l)2(l+0c) + aP} 

Ca 1 )^ ^ 


- C] 



(l+Oj,) + A 

1 



(4 18) 


(4 19) 


(4 20) 


(4 21) 


(4 22) 


(4 23) 


(4 24) 
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1 of 


B = I (a-l) + AKCoK-l) 

r = a-T i to a - 4 C’ - oK) 
3 % 


When K ->■ 0, from egiu1:ionA(4 19)} (4 20) and (4 21 ) the coireBpondlng 
expressions for flow fluX} load capacity and frictional force respect- 
ively are obtained as follows 


0 = 


« 

0+1 


(4 25) 


6 Itn a 

fr . ( . -|£_} (1 + XC ) (4 26) 

(« 1)2 <i2 1 ® 

, 4 to a 

T « {6 _ (1 + xr ) (4 27) 

a+1 o 1 ^ 

The variations of Q, W^} W} (-Tq) and (-t) are shown in 

figures (4 2) to (4 7) respectively It Is noted fromUgrres (4 2), 

( 44 ) and ( 46 ) that (the contribution to liie flux, due to mass 

transfer), decreases while \!IL (the contribution to the load due to 

o 

mass transfer) and (-T^) (the contribution to frictional force at the 
moving plate due to maea transfer) increase as K, the mass transfer 
parameter, increases Erom figures (4 3)j (4 5) add (4 7) it is qbvlous 
that Q decreases and ff, (-T) increase wilh the increase In K 
for fixed Xd 

0 



J OA 



Fie A 2 VAFilATtON Ol( JC 5 OB 




K-— > 




FIG 4 4 VARIATION OF .W# W.HH Ot s ? 2 






K =► 


MG 4 6 VARIATION OF (-Tc ) WITH K FOR tt = 2 2 



FIG 4 7 VARIfTPH OF(-T) WITH ACq FOR VARIOUS \ALUE5 OF K ANDcCs2| 
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Ihe behaviovir of Q, W and (-T) are different with Xo^ 
for each K Sbr example Q ie a conetont vAxen K » 0 and 
decreasee for positive K and Incx eases for negative K vath the 
increase in XO^ [see figursi no (4 3)] But W and (-l) 
increase with XC^ f&r all values of K as shown in flg,ure8(4 5) 
and ( 47 ) 

4 3 HTBROSTATIC BSaRlHG 

Consider the case of a parallel surface hydrostatic bearing 
as explained in figure no (4 8) 2t is aesvsned that the mass is 
being transferred at the plates and also there is chemical reaction 
in the flim Itself, may be due to presence of some detergent In 
such a oaae, tne equation determining the concentration is written 
as follows [eee equation (4 7)3 

u 11= 'C (4 28 ) 

where k' is the rate of first-order chemical reaction and is 
negative if there is a loss in 0, oliierwise positive 

lib obtain simple solution of equation (4 26) we assume that 
the cone entration greidient in the axial direction is large but C 
does not vary very much acrose the film Uhls is possible because 
concentration gredl«it in eoclal direction j?esults from small 
oono entration ohungea over extremely small axial distances, 

SneoJc [1969] With this, taking average of equation (4 28) we have, 







* r 


FIG 4 8 EXTERNALLY PRESSURISED HYDROSTATIC BEARING WITH 
MASS TRANSFER 
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Q II = 2iirC (k'h + k^) (4 29) 

Here q Is tne flow fliuc of the lubricant and Is given by 

h 

Q = / 2irur dz = eons taut (4 30 ) 

o 

and le given by 

-S^h - ll)<. “ ’'i® (-* ) 

Qhis equation represents the mass balance at surfaces SSniB, the 
equation (4 29) may be interpreted as the balance of mass transfer 
in a thin circular fLloi due to convection and due to chemical 
reaction in and around it 

Solving equation (4 29) and using the boundary condition 
C = C^i at r = R , we get 

O' 0 


c 


^(k'h -He,) (r^-S®) 
C 0 


(4 32) 


ashen the average vxscosity of the suspension is written as 
follows » 


U = V„ [ 1 + XC„ Ssp 
o o 


(K> + K ) ( ^ - E^)] 
' R 


(4 33) 


where 
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k.irE 


Now the equation detezmlnlng the pressure In this case Is 


given by 


cl£ ^ 


(4 34) 


where y is given by equation (4 33) Integrating equation (4 34) 
and using the boundary condition P ® at r = and p = 0 at 
r = Ei we get 


5 / r 

irh^ R 


(4 35) 




I 


(4 36) 


She load capacity of the bearing for a constant flux is 


given by 


.-fif / 

h^ R 


yr dr 


(4 37) 


which on using equation (4 33) gives, 


N = 3 (1 - a ) + xa W' , 


(4 38) 



where 



and 

I^iKp {(K'+Ki)(l - R^)}-. 1] (4 39) 

When K' + -»■ (J, we get 

W = 3 (1 -- H^) (l +^C^) (4 40) 

Xhe varlatlonej>f v/' (the contribution to the load capacity 

o 

due to mass treuasfer) e^ W are shown ln£igU7aB(4 3) and (4 10 ) 
respeotiyely It is obvious from these figures that both and W 
increase as E' 4- increases Ihe load capacity also increases 
with the increase in XO for all values of E' + E^os shown in 

0 I 

figure no (4 10) 

Similarly the eocpresslon determining the tor(jue M is given 
by 

2ffft R , 

M = / yr dr 

R 

0 

which on using equation (4 33 ) again ^ gives 

M = ■! (l - R^) + XO^Mp (4 41) 


where 
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and 

V - 5 ^ ‘<''-*1 ’ * Ti^ ' 

(4 42) 

7/heu K' + -i- 0, we lOfAye 

M (1 - (1 +XCq) (4 43) 

Ihe behavlovir of and M are shown in figures (4 1l) and 

(4 12 ) respectively It is seen that both M and ii Increase 

with K • + K. but M also increases with XG^ for all values of 
1 0 

K' + 

4 4 RBSUI/PS AND DISCUSSIOir 

It Is seen ihat the presenoe of additives and impurities in the 

base lubricant is an important factor in detezmining th«^ oharaoterletlos 

of the bearing Xhe presence of mass txansfer process in the 

lubricated system can also modify the behaviour considerably In the 

case of slider bearing it is shown that the load capacity and the 

magnitude of frictional force at the moving surface increase vdth 

XC^; and the mass transfer parameter Xbe flow fluxy however^ 

decreases with K for fixed XC Similar resulte for the load 

o 

oapaoity and frictional torque in the case of hydros tatin bearing are 
also proved* 

finally it may be renarloed tnat if the effect of masa transfer 
(;El or K' + K^) on the surfaces or in 'die film is such that it increases 




-10 -05 0 05 10 

n'*K, — 

V, 

FIG4I1 VARIATION OF VI^ITH K'*K, FOR R- 5 
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the concentration of additives or ImpuritleB In the base lubricant f 
then the load capacity and the frlc tlonal force would increase in 
the above mentioned system 

However^ if the effects of K or K' + are such that it 
decreases the overall concentration of tne additive in the base 
lubricanti then the load capacity and frictional force would decrease in 
comparison to the case of ori^nal concentration Such a situation 
arises when the adaitives get lost on the surfaces or In the film 
because of ohantloal reaction due to the presence of detergent etc 



OriAI!PER - V 


CHARACT:??isa?ics OP iTotr-iTEi’ajoiriAEr Tom iaw f bhicants 
IN AN EXTEBNAliiy ERESSORISED CONICAI STEP BEARING AND 
IN HrUBOSTATIC STtP SEAL 

As pointed out In Chapter the lubrioants whioh do not follow 
Newtonian hypothesis lea linear relationship between shearing 
stress and rate of shear are called non-Newtonian Such non-Newtonian 
lubricants for which shear stress varies as some power of shear-rate 
are callea Ostwald-DeffSale or Bower law models lbs deviation of 
this power index a(8ay) from unity determines the change in Newtonian 
behaviour of tne lubricant Eor n = 1 , the fltdd behaves as 
Newtonian^ for n > 1 it behaves as a dilatant fluid and for n < 1 
it oharaoterlses the pseudoplastla behaviour 

She characteristic e of various bearings with power law fluids 
as lubricants have been studied end the deviatlous from Newtonian 
behaviour have b&en pointed out , Ng and Salbel [1962] , Shukla 
[1963a, 1964b, 19640 ], Hsu and Saibel [1965] , Tanner [1965] , 

Shukla and BreJcash [1969] However, very little interest has been 
shown to study the characteristics of such fluids in hydrostatic 
or hydrodynamic seals, Tanner [1960] , Cheng, Chow and Wllcock 
[1968] 

In this chapter, the uses of non-Newtonian poYfer law fluids 
as lubrloants are investigated in the following systems !Ibe 
effects of stepped film thloKObSe on the ‘various a haraot eristics of 
thfrse systems are polntea out 
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I Externally preseuiieed conical step bearing 

II Hydrostatic step seal 

5 1 BA'IIC lUitTION "OR TH2 PIO ^ OP LAW IBBRICAUTS IH 

A mu PIM 

Consider the flow of a power law fluid in an infinitely long 
thin olearanoe eitb r due to external pre surleation or due to noznal 
motion of the plates In suoh a case the veloeliy profile would be 
as shown in figure no (3 l) 

13ie stress-strain relation characterising the power law fluid 
is given by 


T « 

xy 


m 


* ay' ay 


(5 1) 


wnere m^n are the oonsisl^ncy and flow behaviour indices of the fluid 
respectively Ihe equation governing the velocity of the fluid 
under the assumptiojimof lubrlocttion theory is given by 

(5 2) 

3y dx 

Prom equations (3 1) and (5 2)| we get 


m 


ay 



au, _ t?p 

a7^ 


(5 3) 


Let y e vh be the height where the velocity u of the fluid is 


maximum 


Qhen 
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— = 0 at y = V 1 

a 


(5 4) 


Now in the region vh < y < h, clearly — < 0, and the equation 

— “ ay “ 

governing the velocity of the fluid [from equations (5 l)} (5 2) 
and (5 3)] is given by 


iL. f e . ^ V < y < h 

ay ^ ay - _ 


, V < y 
n dx ”* 


(5 5) 


Integrating equation (3 3) and applying condition (5 4) we have 


I7- (- (y - vh)*/” (5 6) 

3y n dx 

Integrating equation (3 6) a^in and using the no slip oondltion 
u = 0 at y » h, we get the final expression for the velocity of 
the lubricant as follows 


u 


, ~ +1 

i 0-vl-)" 

( ^ 
® dx 1 + i. 

n 


- (y vh) 




vh y £ h (5 7) 


au 


Similarly, in the region 0 ^.y i.vh, we have > n, and the 

3y — 

equation governing the velooily is as lb Hows t 


ay '■ay'' m dx 


C < y < vh 


(5 8) 


Integrating equation (3 8) wi*tih condition (3 4) atvi u a 0 at 
y « 0, we have 
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(5 9) 


1 1 ^ ^ 
u«.( ife" 

m ox 1^1 


^)n _:>: (5 10) 

^ I + i 

1 . y i 

Now at sr a vh the velooltlea are ooatlauous, so from equatloae 
( 57 ) and (5 10 ) we get 


1 1 + 


1 


(- k t 

«* 1 + i 

n 


n 1 + ~ 

(i-v) ^ 


, (- i^)» h - 

" ITT 

n 


i lU 


.V " (5 11) 


where each side of the above equation gives the maxliaum velocity 
of the fluid 

Equation (5 11 ) Inpliee that v » ^ and the flow is aynnsietrioal 
about the middle point of the fll[a*'thxo knees 

Ohe volume flux Q of the fluid is deilned by 


h/2 h 

Q “ t> { / udy + / udy } 

0 h/2 

which on using equations (37) and (b IO) gives 


(5 12) 


Q , 2n^ 1^)11 fij/2) 

^ 2n+1 ^ m dx' 




(5 13) 
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HVhen the flow is due to squeezing, on Integrating the equation of 
continuity ^ ® using the boundary conditions v = - V 

at y « h and v = 0 at y - 0, we have 

a ^ 

V « — / ud/ (5 14) 

From equaxlons (5 12) and (5 14) we hsive the differential equation 
for determining tho flux as follows 


!§* = b V (5 15) 

If there Is no squeezing 1 e V - 0 and the flow is due to 
external pressurlsatloa, Q is a oonstant 

5 2 CASE I EXIT WAIiXf FRBSSURISID CONICAL STEF BE4BIN6 

Consider the case of externally pressurlsea conical bearing 
with stepped film thickness using a power law fluid as lubricant 
The physioal situation Is shown in figure no (5 2 ) 

In this case the equation determining the pressure oan be 
written [from equation (5 13) and putting b >= 2Tr x sin a, and 
Q^= Q (a constant) ]as follows 


dx 


- ^ f (2a+1 )Q ,n/ 2 'j2ntt 1 

4nii Sinn ' * 


1,2 } 


(5 16) 


Where is the flow flux and p^ and pg are the pressures in the 
regions th = h^^k^L <,x ikl) and ^h “ hg, kL l) respectively 




FIG 5 2^£XT£RNAIiLY PR£5St)PlStD CONICAL STEP BEARING 
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The boundary oondltlona are 


P 


1 



Pi 

P2 



at X « k £ 
o 

at X = kl 
s 

at X = L 


(5 17) 


where p is the pressure at the step 

O 

Integrating eqx^tlons (5 16) and uelng the boundary conditions (? 17) 
we have 

p, - T (5 ie) 


and 





(5 19) 


where 


k’-“) + - fc’-“) , 


1 

H«t^ = 
^2 


h-+ h h 

^2 ^2 


(5 20) 


IHirther^ the pres8U3?e at the step and the flow flux are given by 


VF 

Q _ l“ f **1^^ ^ 4o*sina 


(5 21 ) 
(5 22 ) 
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When h-O^le H^l the ooxreepoudlng expression for the flow 
8 

flux is gxyen by 


Q. * £ 




-!■ h 




1 “ ( ~ ) 


2 


4n‘n~eln a 
2n+1 


(5 i3) 


OSie ratio of 7^ with respect to H is shown in figure no* (5 3)i 

^ % 

It oan be seen from this graph that decreases as H or h_ 

Q “ 

Increases and Q ^ Since Is not a function of 

h f It laqplles that Q Increases as H or h_ increases 

S 8 

Ihe load oapaolty of the bearing Is given by 


W » 


2 2 2 

irp.k Ii sin a + / 

'‘o'- 


2irx sin a p^ sin a dx 


+ 


2irklj8inap h oosc-f- 

s B 



2irx sine Pg slnadx 

(5 24) 


Subetl-tuting the -values of Pi>P2 Pg equations (? 18) ^ 

(5 19) and (5 21 ) respeotlvely In equation (5 24)^ we have the final 
expression for the load capaolty^ as follows! 


k^ ain^a + | h^CH-l ) efn« cos 0 

p 

t — sin a 
& 


9 


(5 25) 




89 


where 


^ u _ 2 

c 2 ' 2 L 

TTp^L 


F = {-i k^"“(k^-k^) } + 

c 2 ' o' 


2—} + {1 


4(1-*^) -^!^) • (5.26) 


If k -»■ k y then W is gLven by 
C' c 


3-a\ 2„ 


. (l-k^ ) sin' 

f = ±:a2 — o_J 

o 


+ hgCn-l ) k^ sin 2 a . 


(5.27) 


Ihe aboye result has already been obtained by Shukla [1963 -a] when 
ii = 0, also for a = 90°, i is same as the load capacity of bearing 
with constant film thickness y Shukla S Trakash [l969] . 'lien n = 1 
the load capacity can be written as 


^ . 2 

2 22-3 “ 1 ®^“ “ 

W = k sin a+ -- h.n (H-1 )k Jin k sin a cos a + — 

c O ^ &.] 


(5.28) 


where 


5^ = £n k - Jin ^ 


F = {k^ k " 2 ^^^*^0^^ ^ {(h^-1 ) 

o 


2k !ln k} 
o 


(5.29) 
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Hie -variations of W with H are snown in figures (5*4) to 

c 

• 

(5*9)» It can be observed from these iigures that W increases as 

H or h increases for fixed k and n, and the load capacity of 
s 

the step conical bearing is always greater than the corresponding 
case with unifoim film thicicness [see equation (5.27)1 * 

ilroffl figures (5.4 (5.5 ) and (5.6) it can also be pointed 

out that for fixed H and n, the load capacity can increase or 

decrease with respect to the increase in k and this -variation 

depends upon the -values of H and n. Further from figures (5.7), 

(5.8) and (5.9) it can be seen that for fixed H and k, W 

increases as n decreases for -values of k close to k but it 

0 

decreases for larger values of k. 03aus it may be coxxjludecl that 

for a given n, W attains a maximum -value for some values of H and k. 
c 


5.5 CASE II; fiOH-GONIACIIUG HTDBDSIA'HC SIEP S£AI 


Consider the case of a non-contacting hydixistatic step seal 
as shown in figure no. (5 .IO). a non-Newtonian power law fluid is 
flowing thirough the stepped clearance due to hi^ pressure p^ at 
X = 0. Since the seal width L is small in comparision to the 
outside radius of the seal we can use the infinitely long clearance 
solution for determining the pressure distribution, From equation 
(5.13) we can write the equation determining the presstrre in this 
case as follows: 


dp 

to 


^ f ( 


2 ^2n+1 


h 




5 


(5.30) 
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FIG 5.10 HYC'ROSTATIC NON - CONTACTING STEP SEAL 
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where j = 1 in the region 0 ^ kl and 3=2 in the r^ion 
^ ^ i ^ 3 is flow flxis which is a constant. 

Integrating eq^uation (5.30) and using the boundary conditions 


p. = p. at X = 0 
1 1 

Fl = Pg ai x = kL 

Pg = 0 at X = L 


(5.31) 


we hare the following expressions for the pressures in the two regions 


n = ^ r-p _ r 1 l2a+1 X , - 

^1 


(5.32) 


_ _ / 2 \2n+1 / X\ 

Po -p ^ Vi ' (l ~ _)j kli < X < Ii 


^2 f " h 


(5.33) 


where f = k ( + (l-k) ( 


and h.= h^ + h 
12s 


(5.34) 


!Ehe expression for the flux can be written as 


Q = 


( — } 
^ mif ^ 


1 

n 


2nb 

2n+1 


(5.35) 


When hg= 0, the flow flux in the seal is given by 



2nb 

2n+1 




iio 2 



(5.36) 
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The ratio — is plotted in figtire no. (5.1l) with H for ranoTos 

values of n. It can be seen from there that this ratio decreases 

as H or h increases and C- > Q for H > 1 . again* since Q 
s — o o 

is not a function of h , Q increases as h increases for all n. 

s s 

Hie load caj.cacity is given by 
kL L 

o' kL 

which on using, equations (5.32) and (5.33) gtyes, 


® ^Pi^ 2[k + (1-k) 
Erom eqmtion (5.37) we have 


(5.37) 





k(l - k) 

2[k + (l-k) 


> 0 


(5.38) 


Since H > 1 and n > 0, IT increases as H or n increases. 

Thus we observe that f increases as H or h increases for 

s s 

fixed k and n. Further, it can also be concluded that W„ 

' s 

increases as n increases for all H > 1 and 0 < k < 1 . ’See 

variations of W with H and a for different values of k are 
s 

shown in figures (5*15) to (5*17). Again, tTcm equation (5.37)» 
we obtain 




VARIATION OF QoJ'Q' WITH H IN ‘fHE'CASE OF A STEPSEAL 
FOR k =0._5. : ' , , . 
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aWg 

It is seen that sigi of -r^ depends upon H and n for fixed k 

and is zero for some value of k for given n and H. Hence Wg 

can increase or decrease with increase in k depending upon the 

values of H and n. Ihs variations of 7 with k can be seen 

s 

from figures (5.12) to (5.14). 

Ifh = 0(H = 1,k = o), from, equation (5.37) we have 
s 


Wg = 1/2 (5.39) 

It IS noted from here Ihat the load capacity in case of a seal 
with uniform gap does not depend on the flow behaviour indices 
of the lubricant. 

5.4 COFCniSIONS 

In this chapter, we have studied the effects of stepped film 
thici-cness on the various characteristics of the bearings and seals 
using power law fluids as lubricants. 

Ihe following conclusions may be drawn from the analysic: 

(1) In the case of esrsemally pressurised conical bearing the 
flow flux increases due to increase in the step hei^t. Similar 
situation occurs in non-contacting hydrO-static step seals. 

(2) Ihe load capacity of a conical step bearing is always 
greater than that of a corresponding bearing with uniform film 


thickness 
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(3) In the case of a conical bearing, fbr given H and k, the 
load capacity increases as n decreases for smaller values of k 
but the situation is reversed for larger values of k. itirther, 
for given H and n, the load capacity may increase or decrease 
as k increases and this variation depends upon the choice of H 
and n. 

(4) In the case of hydrostatic step seals, the load capacity 
increases as H or n increases for all values of k. !Ihe load 
capacity may decrease or increase as k increases, depending upon 
the values of H and n. However, after certain transitional 
value of H depending upon n, the load capacity may increase with 
increase in k. 



GHAIT-ffi-TI 


fHSRMAl EFi-JG^S IN EXTSEaSTAIir IHESSD'RISSD BSABING WITH 
PO’SER lAW IHBRICANTS 

la hi^ temperature and pressure applicatioas the flow behaviour 
characteristics of the lubricant (Newtonian or non-Newtonian) no longer 
remain constant, but they vary with tanperature and pressure, Rotem 
and Shinnar [1962] , Tuzlan [1965] , I)yer [1969] . Several 
investigarions have been conducted in recent years to study the 
characteristics of bearing by considering these vacations in the case 
of Newtonian lubricants » Ting and Mayer [1971] » Rodkiewicz and 
Anwar [1971] # Donaldson [1971] t Go^d- [1967» 1971] • In particxilar, the 
effect of temperature in externally pressurised bearings has been 
studied by Donaldson [1971 ] , Ting and M^er [1971 ], using exponential 
relation between viscosily and temperature.- 

In this chapter, we study the characteristics of actornally 
pressurised bearings using non-Newtonian power law lubricants vsfaere 
consisterKsy of the lubricant varies exponentially wilh toaperature 
and pressure. 

6.1 BASIC EQUATIONS 

As before, let us consider the flow of a non-Newtonian power 
law liabricant in an infinitely long thin clearance due to external 
pressurization. llie physical situation ard the co-ordinate system 
are shown in figure no. (6,1 ). 



y 


/ Z ^4-^l.Z..Z.Z „ Z Z ..Z.Z..Z.. Z -^ " 

^ t 

h/2 

t 




0 


•^u 


->x, 


o' 


/T"7 7 ■/■■/~7-7-7-rn A ^yy.y-^.y-~-r-y~V- 


/I 






1 -/ 2 - 


FIG.61 EXTERNALLY PRESSURISED RECTILINEAR PLATES WITH 
POWER LAW LUBRICANTS 



FIG 6 2 EXTERNALLY PRESSURISED CIRCULAR PLATES WITH 
POWER LAW LUBRICANTS 
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Ihe equations of momentum and energy wnich govern the flow 
of a power law lubricant are given by 


C 



JL 

ay 


ay 


aui 

ay 


( 6 . 1 ) 


Cy ^ 


dT _ 


dx 


= k 


ay'^ 




n-1 


aui 3u 

ay^ ay 


( 6 . 2 ) 


where m, the consistency index, varies exponentially with temperature 
and pressure as: 


m = m 

o 


e 


ap-3T 


(6.3) 


Since the physical configuration is symmetric, we consider the 

region 0 5_y 1 . fUrtner investigation. In the region 

0 < y ^ht/2, clearly <0 and equation (6.1 ) can be written ais 

ay “ 


3 ^-,11 _ . dn 

ay *■” ay^ " ~ n 


(6.4) 


Integrating equation (6.4) and using the boundary condition 


we get 


■p = 0 at y = 0 

ay 

u = 0 at y = h/2 , 


an 

ay 


_ ( _ 1 ^ V 

^ m dx ^ 


n 


1 

n 


(6.5) 


( 6 . 6 ) 



112 


u 


'' m 




1 + 



- y 

a 


31ie TOliaae flux of the lubricant is defined by 


(6.7) 


0 h/2 

'^ = b {/ , udy + f udy} 

-hr. 0 

which on accoimt of symmetry reduces to 


li/2 

r = 2b / tidy . (6.8) 

0 

'^ing use of ecLuation (6.7) in equation (6.8^ we obtain the final 
expression for flow flux as 




2nb 

2n+1 



(6.9) 


It can be easily seen firom the eqttation of continuity that 
the flow flux remains constant and does nof depend on x. In order 
to consider the effect of tauperature we write the equation (6.2) 
in the following form: 


pc k n (- , 0<y<h/2 . (6.10) 

^ rtx oy — 

AreiBging the equation (6.1 o) from y = 0 to y = h/2 aM making use 
of the symmetry we get 
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r dT . . Du-r+l . 


( 6 . 11 ) 


3 OP “i" ii /2 

where ^ 1 _ J 3/2 nteglectfcd. 


Irom equations (6.6), (6.9) aod (6.11 ), we finally obtain. 


pc 


dT _ <^P 

V dx dx 


(6.12) 


!I3ae boundary conditioas for T and p are 


1 = 0 at X = X 


p = p. at X = x„ 
i c 


(6.13) 


Integrating eqiiation (6. 12) and using the conditions ( 6 . 15 ) 


we get 


T = - 7 (p - Pi) 


( 6 . 14 ) 


where y = pc. 


V 


ilrom equations ( 6 . 3 ) and ( 6 . 14 ) we get the final expression 
for conffiLsteney mg tax sis 


m = Wq e 


1 XT 

Y JI ^ 


(6.15) 


/ 8 V 

■tiiere x=(a+ ~/Pi is a diaiensionleBS parameter. 
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!Ehus the depend ence of consistency on temparatiire and pressure 
ultinately reduces to pressure dependence in this case. 

6.2 ESaSEITiiLn’ PRESSUEI5SD INPINIOSIf iXTJKDED aSCiPiinGULxiH OHRUSS 
PlaTES 

In this case, from equations (6.9 ) and (6.15) the final expression 
for determining the pressure of the lubricant can be written as; 


2nb 

2n+l 


8Pi 




^ 1/n 

dx"^ 



a constant . 


(6.16) 


Ihe boundary conditions for p are 


p=p. at x = x , 

^ ^1 o’ 

p = 0 at X = • 

Integrating equation (6-16) and apply’ing (6.17), the final 
expression for p is given by 


(6.17) 


p = - Jin [ 1 - ( i:^) (l - e"^)3 , (6.18) 

* 1-ac^ 

X 

, “ X - o 

where x = -r » x = — . 

L o li 

Erom equations (6,16) and (6.18), the expression for the flow flux 
can now be written as 


(6.19) 
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Q 


_ , 2v.il - e”^) e ^ , 2 + — 

2db |- -1 n 


2n+1 


L XCl-2x^) 


(6.19) 


!33ae flow flux for constant consistency index is obtained from 
equation (6.19) by making X j g 2 as follows: 


0 


_ 2nb ^i -.n ( ^ ^ 

2n+1 ‘■m 1(1-2^ )■' ^ 2 ^ 


1 1 
2P. ^ . b 2 +4 


(6.20) 


Erom equations (6.19) and (6.20) we get 


^ = E- 

r 


^ (1 - e-hj H 


( 6 . 21 ) 


!I3ie locid capacity of the bearing is cLefined as 


Xjj L/2 

W_ = 2b / p.dx + 2b / pdx 

X A X 


which on tising equation (6.18), gives 


( 6 . 22 ) 


n_ 


Wn = 

P b -n . ! 


1 - (1 + X)e"^ - 2x^(1 - X - e"^) 
XCl - e-^) 


( 6 . 23 ) 


It is interesting to note frcsa equati£>ns(6.19) and (6.23) that for 
prescribed applied pressure the flow flux depends upon the flow 
behaviour characteristics of the lubricant while the load capacity 
is independent of these indices. 
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6 3 SXTiJiNALEr HlBSbimiSSlD ClRCULiia 3i!.iiRnrG 

In this Bsction^we consider the case of externally presaunsed hearing) 
the physical coaflg;ur>atlon of which Is lllustr-ted in figure no (6 2) 

In this case) from equations (6 ^ ) and (o 15} hy putting b = 2icrf 
the final expression for detexalolng the pressure of the lubricant 
Is given by 


(- i_ e Y e 
2n+l Pq 




n j ~ "> + i 


?s 0 = a constant (6 24) 


The boundazy conditions for p are 


p = Pi 

at 

r = r^ 

p « 0 

at 

11 

u 


(6 25) 


Integrating equation (6 24) and using the condition (6 25) 
we get the expression for p as 


p.3 ,n [1 -(1 -e-M 

^ 1-^0 


1 


(6 26) 


where ;=|, = f 

Erom equations (6 24) and (6 26) the expression for the 
flow flux Q is ^yen by 

^ “ot’’* ” - *o'"5 


(6 Z7) 
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Here also the flux corresponding to the constant consistency 
Index Is obtained by letting 6 and X *> 0 , as follows 

(6 28) 


(6 29) 

9!he load oapaoii^ of the bearing Is defined as 


1 


1 


I f 


2n+1 




From equations (6 2?) and (b 28) we have 




^ r - 

(T" ^ ® 

'• o 


1 e-\n 


U = w r^ + 2ir / rr dr (6 30) 

which on using equation (6 25) gives 

y - — ! — — s _ ^ ^ dr (6 51 ) 

imi r^ dr 

1 o 

Substituting for ^ from equation (o 26 ) ^ the final expression 
for load capacity is given by 


VI s 


(1 e^^l (I n) 1 


-2-*n 

T 


^o Cl-rJ^'D-a e-^)(l-r^’''3 


dr 


(6 32 ) 


ihe expression for the load capacity for Kewtonian lubricant is 
obtained from equation (6 32) by letting n -> 1, which gives 
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t; 


1 e 


X 1 


L 


T dr 


^ ^ fl An r An 


(6 33) 


In the ease of oonstant consistenoy Indox the equations 
(u 32) and (6 33) slmpliiy to (by loakiiig X 0 ) 


o 

^1® TT » n - 1 

^ 2An(i ) 


(6 34) 


(6 35) 


It l8 easily seen that equation (6 34) is the same as equation (3 27) 

with h B 0| while equation (6 35) is the well known result for 
s 

hydrostatio thrust bearing y finkus [i960 1 

She expressions for load oapeoliy for large X can be written 
from equations (6 32) and (6 33) ast 



"i 



L 


1 

I 


(1-g 

- l-n - 1-n 
r — r 
o 


dr , 



n / 1 


(6 36) 


It is easily ^een from these equations that for large X the load 
capacity deoreases as X Inoreases for all n 



119 


6 4 Discussion sm bssui/es 

It Is seen from equations (6 21 ) and(o 29) that the expressions 

for Q/Q^ In both the oases of rectxlinear and circular plates are thb 

same She variations of q/Q^ are shown in figures (6 3) and (6 4) 

l*P, 

for different n It is noted that Q/^ increases as 

increases for fixed It can be further seen from these figures 

that the ratio Q/%, can increase or decrease as n increases and 

this change depends upon the choice of ap^ and 

1 Y 

She variations of Vp In the oase of rectilinear externally 
pressurised thrust bearing are shown In figure no (6 6) for two 
value e of reoess It is noted from Hie expression for Wp from 
equation (6 23) that the load capacity is independent of the flow 
behaviour indices of the lubricant :Bhrtber,lt may be pointed out 
that Wp decreases as \ increases llhe variations of the 
dimensionless load capacity In the case of externally pressurised 
circular plates are manifested In figure no (6 5) It may be 
mentioned that V decreases as X or n increases 

Dsr putting n = 1 and x^* 0 or r^» 0, we get the results 
obtained by Donaldson ri97l] Burthert when X 0 ^ we get the 

corresponding result of non-Newtonian power law lubricant without 
thermal effects [see equation (6 34)] » [Sbulcla 1963**a] 





1 



FIG 6 5 VARIATION OF Wp WITH A IN CASE OF EXTERNALLY 

PRESSURISED RECTILINEAR KATES FOR DIFFERENT Xq 



FIG 6 6 VARIATION OF W WITH 
PRESSURISED CIRCUL 



OHArOJiB - VII 


IXIJRViililiSr IRiiS&UBIbJiJ} K>BOUS !|!HRUST BEAHINO WIXH 
SOW^ lAV UJBRICMiTS 


Sie oharacteilatlos of externally pressurLsed poroue hearinge 
uelng Newtonian fluide as lubricants have been studied recently due 
to tbelr much wider applications i Iforif et al [1965] t Hslng [1971] 

As most of the lubricants are polymer solutions ^ tlie question natvirally 
arises as to how the oharactexlstlcs of bearings change when such 
rheological substances known as non^Newtonlan fluids are used as lubricants 
Oils and other related questions have been answered by many authors who 
have studied the oharaoterlstloa of various non’-poro\xs bearings by 
considering Iow«r law models of the lubricants ^ Ng and Salbel [1962] t 
Tanner [1963] » Sbukla [1963a] , Hsu and Salbel [1963] t Shukla and 
Irekash [1969] In particular , the oharaot eristics of power law 
fluid as lubricant have been studied in externally pressurised conical » 
Shokla [19b3a]and olrouiar thrust bearings^ Sbukla and Trakash [1969] $ 
without taking pecoslty into aooount 

In this ohupter^we study the ebaraottoxistios of externally 
pressurised porous thrust bx^rlng using power law fluids as lubricants 
Here we have assvaaed both the plates to be porous to exploit symmetry 
In oalculations, howeveriln aotual applications only one of the plates 
needs to be porous and the results presented here would be equally 
applicable in ihls case iSie analysis la based on the assumption 
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that the poroua natrlx conalsts of a syetam of capillaries of 
very small radii and thus restricting the flow of the lubricant 
throu^ the matrix In only one dlrectloni Hslng [1971] As 
the exacc analytical solutxon of the problem is not possible ^ 
both approximate and numerical solutions are presented She 
physical situation of the system oonsldered is shown in figure 
no (7 1) 

7 1 BASIC EQUaXIOBS 

Ihe basic equation governing the flow of ihe lubricant 
between the bearing surfaces » under usual lubxioaplon assumptions 
( as in 9 5 1 ) is given by 





where 


^xy 


.au.n-l 

” 'lyl 


au 

3 / 


(7 1) 


SlnoS) in Hae region H ^y ^ H + ^ i from equation (7»1 )» 

we have 


K A" , i ^ 

3y &y in dx 

Integrating equation (7 2) with boundary conditions = 0 

ay 

at yaH+^<ussO at y = H, we get 


(7 2 ) 



|24 



FIC 71 FLOW OF A F OWE(7 LAW FLUID IN A THIN fLEARANFF 
BETWEEN TWO EXTERNALLY PRCSSUM1£D TaRALh-L 
POROU ■) I I ATES 
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u 


1 

f _i iE 

' m dx ' 






(7 3) 


Slmllarlyi in the region H+^^y^H + h, — <0 and the equation 
governing the velocity is as follows 


3y m dx 


(7 4) 


Integraclng equation (7 4) and applying the followxng conditions 


1^=0 aty = H + V2 
u^O at y = H+ h, 

we get 


u = ( 


lifi \1/n 
m dx ^ 


[(| )’*“-{y-(H +!)}’*“] 

H-l 

n 


(7 5) 


Now^ on integrating the equation of continuity 


9u 

ax 


+ 



(7 6 ) 


(-v)l 


H+h_ 

H 


a 


,H-»h 


/ udy} 
H4h/2 


we get, 
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which oa using equations (7 3) and (7 3) gives 


2 + ^ 

. h V n n / I iE OA ^ 

' 2 ' 2n+1 dx ' m dx ' ’ 


(7 7) 


where v' is the average velocity with which the lubricant Is 

Cl 

injected into the £L1 di 

7 2 MOSIHBD lOM OP SABCY'S LiiW 

Let us consider that the porous matrix consists of a system 
of capillaries whose axes are directed tovrards the film !lhe 
motion of the lubricant In a typical capillary In the lower plate 
is governed by [see figure no (7^2)] 


- i 1— (r T* ) 

dy r 3r ^ ry' 


(7 8) 


where 


T' 

ry 


Dv^n 


3v* 


) since ■! "* < 0 

gy /jaiW'.Wgj. „ 


Integrating equation (7 8) £u3d using the boundary conditions 


3v' 

ar 


= 0 


V*ss 0 


, r =0 
, r - E 


(7 9) 


we get 


yt 



n 


(7 10) 




FI6 7 2 aOW OF A fWER LAW LUBRICANT IN A THIN POROUS 
CAPILLARY 
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The flux miiy be deUxma as 
P 

O' a / 2ir rv' dr , 
0 


which} on using squatxon (7 10)} gives 


^ , .2m (7 11) 

dy 3+i. 

nir R " 

where Q' is a constant as can be seen from the eqviation of 
continuity 

Integrating equation (7 1l) and using the conditions 

p* = p at y «=: 0 and p' = p at y = H, we get 
0 


Q' = 


nir 

3n+1 


‘i;r> ® 


(7 12 ) 


^e average velocity v' of the lubricant can then be given 


v' - ^ 

^ we 2 3n + 1 


n-p i 

( 

^ 2nH ' * 


(7 13) 


where is the average radius of the system of capillaries 3his 

a 

equation m£^ be interpreted as a modified fbrm of ]3aroylB law Ibr 
power law fluids 

7 3 HBlDEBMlSAnOS 07 lOiL 

Srom equations (77) and (7 13) we obtain the basic equation 
governing the pressure in Hie film as follows: 
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I b I I 

r « JL ^ - 2n+1 ( 2 \ / ^ -Jn 

dx m dx ' 3n+1 ' h ' "a ' gmH •' 


(7 14 ) 


])Ub to symmetry I the boundary conditions for the pressure are 


dx 


0 , 



X = 0 


p = 0 , 



(7 15) 


fdiere p 
ni 


Is the nkaxlmum pressure 


Intzoduoing liie variable C * T“ 

It 


£ s p - p in equation 
& 


(7 14) we get 


1- 

d5 W 


h ^ 


n 


(7 16 ) 


where Kg 


4H / 2n+1 V / ”a^ J'*' n 
h ^ 3n+1 ' ' 2Hh 


(7 17) 


18 a dimensionless parameter depending upon the oharacteristlos 

of the porous matrix and p is the externally applied pressure 

© 

Integrating eqmtion (7 16) and using the conditions 


45 


0 and J = at 5=0 

HI 


we get 


45 ^ 




(7 18 ) 
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and 


P 

/ 


dP 

1+i uj. 

fP ^ - P 
^ " m 


n 

(n 


€ 


(7 19) 


Here = p^- and can be deteanalned by the followoiig equation 
which 3 .S obtained by using the boundary coodltlon P = p at £ - 1 
In equation (7 19) 


^ rr' ' 


She load capacity W is given by 

L/2 i 1 

W « 2b / n dx » b L / p dC = UL / (p P) d5 , 
^ 0 0 


which on using equation (7 18) gives 


W * 1 - 


2ft o’B’ 


/ 

1 


z dz 


(! - 1 ) 


n 

n+1 


where W 


V 


bip, » 



z 

0 




(7 21) 



2 

2n |n .. n+1 
3n+1 



and a’B' 
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Hablly, the equation (7 19) for datexmlniog with the above 
traneformation could be written ae follows 


/° -i* — — = a'S' (722) 

1 1+i JL. 

(z " 

7 4 DIS3USSI0N and RESULTS 

Ib determine Ifae load oapacity of the bearing, first the 

values of Z. for different sets of values of a’ and s' are 
o 

numerically computed from equation (7 22)and these values are then 

substituted to determine V numerically from equation (7 21 ) The 

H 

variations of W with 0' for different values of ^ and n are 
shorn in figures (7 3) and (7 4) Erom figures no (7 3)> (7 4) 
and (7 5 ) it can be seen that the load capacity increases as n 

JJ 

increases for fixed 0' and ^ EUrther, from these figures, 

it is clear that for fixed n, the load oapacity inoreaees as 0' 

H 

and ^ increase i e as the ohorao texistics of the porous matrix 
increass Xhese results oan be seen from the Jbllowlng approximate 
calculations also 


Let us assume that the externally applied pressure is large 
1 e P- » p With this assumption, equation (7 14) becomes 

w 


d / 1 dp\l/n 2a+1 f 2 \n 

dx m dx^ " 3n+1 ^ h ' "a ^ 2mH ' 


(7 23) 


following the earlier procedure, the load oapaolty ui^er this 






i 
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approximation can be written as follows 


w - f iS / gp+1 f- JL 

® ^ h ^ ^ 3STT ^ ^ 

- 4H 

It is obvious frcM here that W Increases as p' or increases 

Furthte?! it can be seen that W also Increases as n increases for ^ 
4H 

and greater than unity 

rlhen n - 1 , i e the case when the lubricant is Newtonian) 
analytical expression fcr V is obtained as 

tanh $' a. 

u s 1 - ^ 

S' Oj 


where 


,3Ha/2 


which is the same as obtained by Hslng [1971] when inertia effects 
are neglected in his case 



CHAMJR - VIII 


SQUEEZE EIIMS USING lOVIR JM lUBRIOi^TS 

Id. the previous three chapters we have studied the behaviour 
of power law lubricants In externally pressurised bearings and the 
effects of step in the film thickness^ wriatioo of consistency 
index with respect to temperature) and porosity of the bearing surfaces 
have been investigated 

In this cfaapxer)We stitdy the characteristics of squ eze film 
bearings using power law lubricants by tddclng Intertla effects and 
the linear variation of conslsxency with pressure Into account !l!he 
effects of step In the geometrical film thickness are also studied 

In recent years ) several Investigators have studlea the e foots 
of inertia in various bearings by oonsiderlng Newtonian fluid as 
lubricant) Osterle and Saibel [1955] ) Eowson [1961] , Milne [196?] > 
Snyder [1965] ) Kumna [1967] f Xing and Mayer [1971] Xhe effects 
of viscosity variation with temperaturs and pressure have also been 
studied In the case of equeeze film bearings Gould [1967] !I!hou^ 
the charaoterlstlcs of squeeze film and other bearings using power 
law lubricants have been studied Ng and Saible [l962] ) Shukla [1964^] 
Sbukla and Erakash [1969 ] ) effects of inertia have not been taken 
into oonslderation In view of tblS) we study the following effects 
in the squeeze film bearings) 
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(1) j£ffeot8 o± Inertia and pressure 

( 2 ) ilffecte of step la the flM thickness^ 

!Ibe cases of Infinitely extended rectangular and circular 
parallel plates are considered 

8 1 8QUi£BZB FIIMS BEirViri&M«r inO INMNlZb. 

Let us consider the flow of a power law lubricant between 
two Infinitely long rectangular plates, the physical configuiatlon 
of which Is shown In figure no (8 1 ) We assume that the two 
plates are approaching each other with relative velocity V and 
the lubricant Is being squeezed out due to this relative motion 
Making use of the symmetry of tiie physical situation, the basic 
equations governing the flow of the lubricant can be written as 


3 X 


3y 



3> 


{-B(- |3") 


(8 1 ) 


3u 

8x 


+ 



lH.<rin''<y< h/2 

* ay ~ ^ 


(8 2 ) 


Here the oonelstency Index m le assumed to vary with 
pressure linearly as fbllows 

m a (1 + op) (8 3 ) 

where and a are constants Uhls relation is compatible with 



• • 



FIG 8 2 SQUEEZE FILM BETWEEN TWO CIRCULAR KATES 
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the exponential yarlatlon for email values of a in which we are 
Interested Since the fluid film is very thin, on taking the 
average of inertia terms equation (8 1 ) can be approxiiik»ted as 
follows SlhkuB [1960] 


1- 

ay 


av^ 




2o 


As p is a function of x, the rl^t hand side of equation (6 3) 
can be written as a function of x and thus we have 


1 - 

ay 


C 


ay' 


* £(x) , 


(8 5 ) 


where f(x) 


& 


2p 


m^Cl+a’^) h BipCl+OT') 0 


h/2 

/ (u + V dy 
* ^ ax yy* * 


She boundary conditions for u are 


u - 0 at y = b/2 ^ 

•|^ = 0 at y * 0 

ay 


(8 6 ) 


Solving equation (6 3) and tislng the conditions (8 6) we obtain the 
expression fbr u as 


u 


1 




u 1+ ~ 

(I) “ 


" y 


1 + - 
a 


1+1 

n 


(6 7) 



Substituting Ibr from equation (8 7 ) In equation (8 2) and 
using the boundazy conditions 


< 

n 

o 

at 

o 

11 

V 

^ “ “2 

at 

h 

y = 2 ' 


we get 


V = 



1 




} 


and 



Integration of equation (8 10) gives 


(8 8 ) 


(8 9) 


(8 10 ) 


(f (8 11) 

when o is a constant to be determined 

FixMtt equations (8 2), (8 5), (8 7) and (8 11 ) we get 


- op) (^V)"(lr«)“ cl )“^' 


2n 


+ 4p 


2n+1 

3tH2 



(8 12 ) 


fhe bovtndary conditions for p are 
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P = 0 at X = Lt I 


• = 0 at X = 0 


(8 13 ) 


Using the second condition of (8 13) and keeping in view the case 
n » 1, we get c ~ 0 Then equation (8 12) givesi 


Q 

dn . \ / 2n+1 „\n , 2 \2n+1 n. 4p(2n+l) V 

. J = m^<H-ap) ( — T) (^ ) ^ 


(8 14 ) 


Integrating equation (S 14)fand now using the first condition of 
(8 13 )f we have the final expression fbr p as 


2 

^1 *^2^2 _2, . 1** ,-n+1 n+K2 

„ .-I t I* “ X /+ A 1. 1< X / 

“ ^ ^ 2(n+1 r 


a. a a 
12 


(8 15 ) 


Where a and its hl^er powers have been assumed negligible in 
obtaining the above solution, and 


‘1 = “o ( 


2n+1 xn / 2 x2n+1 , 
2n h ' 


2tt+l\ 




The load capacity W for W3^th b la defined by 


(8 16) 


W a 2b 


(8 17 ) 
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which OQ using equation (8 15) gives 


ft * 2b [ 


1 


n+2 


d L 
2 


2 ^2n+3 

a. aJ, 


d-d al 


n+4 


n+2 


(n-t2) (n+3) (S;!) '(nH)^ 


Sie corresponding case wl1h no pressvtre variation oan be derived 
by putting a s 0 in equation (8 18) as follows 


(e,5) 

h 

Since this equation is too oomplicated to sol-ve for V we determine 
the time for squeezing by neglecting inertia but taking pressure 
variation into acoovuat Q!hen equation (8 18) oan be slnQ)llfled to 


where 


W = K^o©^ + KgB^ , 

B K - 2b(n-H 

1 n+1 ^ 1 (jj+g) (2n+3) 


and 


_ ai)(n-n )l°'^ 
2 n+2 


(8 20 ) 


(8 21 ) 


Solving the quadratic equation (d 20) for and retaining terms 
of the order of a ^ we get 


B 


1 




(8 22 ) 


How putting V « -|t dnd inliegrating equation (8 22) we get the 

d V 

expression for the elapsed time required for the film thicknees 



145 


to be reduced from to for a given lomd V^y as 


. 2 ( 20+1 ) 
n+1 


4m„b 

/ o 

wTo^) 


i 



1 




•}{1 + 


(n+2)Wa 1 

2n(2n+5)bL 

(8 25) 


When a = Oy we get the r.,sult <4S obtcdned by Shiikla [I964b] 

It can be seen from equation (8 18) that the effect of Inertia is 
to Increase the load capacxty of the bearing It can also be seen 
that the load capacity VT Incireasee as a locreasee When a - Oy 

since ) is positive y It Is obvious from equation (8 19) that 

the contribution to the load oapacity due to inertia Increases as n 
Increases 


Erom equation (8 25) It is easily seen that the time of 
squeezing also Insreasea as a increases but the contribution due 
to a decreases as n InoreaseBy since ^li(2n^7^ negative 

8 2 SQU JZE EEIMS BSTWEEJT a?W> CIRCUIaH PlA3?l!.S 


Consider two circular flat plates of radii R approaobing 
each other with relative velocity T !Qie clearance space between 
the places is filled with non-Newtonlon power law lubricant which 
is being radially dlsplaoed outward by the relative motion of the 
plates JSie physical oonfiguratlon is Illustrated In figure no (8 2) 


Taking into account the inertia eifeotSy the basic equations 


governing Hhe flow pf the lubricant In the upper half of the film are 
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iL 

az 




dx 8r a* 


ifS'l 0 in 0 < 2 < h/2 


(8 24) 


Eere also we assume That th^ consistency Index m varies linearly 
with pressure as given by equation (8 3) 4galn averaging Inertia 
Terms across the film thickness ^ from equations (8 3) and (8 24 )| 
we have 


^ (- e ,Kr) , 

3z 32*' ^ ’ 


(8 26) 


where (ftCr) 


dn 


2p h/2 

/ (i^ly+vl^dz (8 27) 

ra^Cl+op) n^hCl+on) 0 3^ 3 


Solving equation (d 27) with the boundary condltionss 


u B 0 at a s h/2 t 


ds 


0 at a = 0 


(8 28 ) 


we have 


u « 4 


h ^'*‘n ‘*'*‘n 

(|) “ 


(a 29) 
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Now subatltutlng u from equation (8 29) in eqmtlon (8 29) and. 
integrating with following boundaxy cooditlonet 


V = 0 at z K 0 


v = -|at a- 1 


(8 30 ) 


we have 


and 


ul 


2+1 
n 


, (3,1) 

14*— 2 +-^ 

a 


1 p. i 

1 3_(r4“ ) , )T ( 2 ) “ 

r dr ^ ' 2n ^ h ^ 


(8 32 ) 


Integrating again eqmtloa (6 32) we get, 


^ _ r (2n+1 )v,n ,2 , 2n+1 ^ a 


(8 33) 


where A ie a conatant be detemlued 


Broni equatxona (S 25), (8 28) and (8 30) we cam wxite 


(- I? * V |H) dz . /‘'\2 u|H* ^) dz 


ar r 


(8 34) 


whiob on using equations (8 26), (8 27), (8 29) and (8 34 ) gives 


dp sr(2n*l)V,nA2n»l^^ ^ ^ £. 2m-l ^ ^ 


C8 35) 
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Xhe boviudaiy condltiona p are 


•d 

fl 

o 

at r = R 

|£ = o 
dr 

at r = 0 


(a 36) 


She second condltaon of (8 36 ) when \ised in (6 33) Implies that 
A = 0 and we have) 


dp , (2n+l)v ,a / 2 s2n+1 a 3 P 2n+1 

-^ = m^(l+ap)[l_^] (j) ^+T3S^7 


(8 37) 


Integrating equation (8 37 )» using the first condition of equation 
(8 36 ) and retaining terms upto a ^ ^ final eicpression 

for p as 


^ r"’) + ^ (H=- r") SJLj 


n+1 


2(n+1 ) 


^3*^4 « 


2(n+l) (a+3) 


{(n+3)R^(R“‘"^- r“+' )-(n+l)(R“+5. 

(a 38 ) 


a.hor« r. On, ,(2a+l)Y U . 1 v2n+1 

vthere a = 2m^ ] ( ^ ) 


(8 39) 


and 


a 


ale, JS^ 

’4 ” 2 3n+1 jj2 


Obe instanteneous load oapaoity W for a given squeeze velooity 
is given by 
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p 

’ * / 2irr tj d-** ^ 
0 


which on using equation (8 38) gives 


I » 2ir [. 


3 a4p^ aj a 

+ + .. — — . + 

2 (n+3) » 






] 


^(n+2)Cn+3) 2Cn+3)Cn+S) 


(8 40) 


(S 41) 


!Ehe load capacity for constant consistency Index is obtained by 
putting a s 0 in equation (8 41 ) as follows 


V « ^ ^211*1 y^A^1^2n+l pn+3 ^ ^ 2nn £_ ^4 


n+3 


n 


8 3n+2 ^2 


(8 42) 


liquation (8 41 ) is the general non-linear equation fbr detexminlng 
V and henoe the squeezing time Sowever f to see the effect of a 
on the squeezing tlnB| we neglect the inertia tezms in equation 
(8 41 ) and we then have f 


+ K^Bg - W *= 0 


(8 43) 


where 


and 


B. 


n+l ’ 


E at •!£• 

*■3 2 


(a4l )^ j^2n+4 ^ 

(n+2)(n+3) 


K 4 =» 


i>(n+l) ^+3 

n+3 


(8 44) 
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Solving equation (S 45) for and retaining teias upto a onlyi 
we get 

Bg =/~ (1 - ^ o ) (8 45) 

ami then as before > the squeezing time oan be written as follows 


t 


a 1 , 3 

2n»l p ^ % , n p n, 1 

Cn+3)W^ 1 + i 

ho ” 


*1 

n 


on 


W n+3 
2n Ti(n+'^) p2 


•] 


(8 46) 


We got tho results already obxainea by Sbukla [1964b] when a » 0 
Aa in the pravinua section) it can be aeen from equation (8 41 ) that 
the effect of inertia is to inorease the load capaci'^ It is further 
noticed that the load capacity increases os o increases In. the 
absence of inertia term) thu squeezing time increases as o increases 
and this xncrease decreases as n Increases 


8 3 EIFECTS OF ST&f IN Flltf TH2DKNBSS 

In this section the effects of step in the film thickness of 
squeeze film bearing are studied She cases of two infinite plates 
and two oirculsr j^tes are coneidexedt 

A SQUEEZE HIM BEOTWEEIN fMO ^TEBiED fARATIEL PLATES 

Consider that a power law fluid is being squeezed out in a 
stepped olearance between two infinite plates due to nozmal motion 
of the upper plate as shown In figure ho (^3) 




FIG a ^ SQUEEZE FILM BETWEEN STEPPED PARALLEL CIRCULAR 
PLATES 
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In this ca8« the fluid flow is synmetrlcal with respect to 
the line x s o aol the film thickness in the two regions are given 
by 

h = h , 0 < X < kL 
h = hg , kL < X £ L 

If Qod P^)P2 volume fluxes and pressures 

respectively in the Isvd regions ^ then from equation (JS 15) we have} 



bv 


where j = 1 for 0 £ x ^ kL 
and j = 2 for kL X _< L 


Integrating (8 48) we get 


(8 48) 


Q. = bVx + A 
J 3 

Since Qg at x = kL} we get 

Vj^=bVx+A> t'^x£kL 
(^sb\x + A kL^x^L, 

vdiere A^= a is a oonertant Ihis constant is zero as 
at X = 0 

ShuS} finally cun be written as 


(8 49) 


(8 50) 


0 


Qj = b 7x (8 51) 

Pcom eqtjstions (5 13) saad (8 51 )$ tiie equation detezmining the 
pressure is given by 
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/ 211+1 „xn . 2 \2n+1 a ^ _ 

to =- ( h ) == . 3 = 1.2 

9 

lategvatiQg this i^nd using the boundary conditions 


(8 52 ) 


P., = Pg at X = kL 
pg = 0 at X = Ir 


we get 


Q / 2n+1 2 \2n+1 r /». _ \n+1 n+1 \ . / .2 \2ii+1«n+1 / . i o+1 \i 


° / 2a+1 \n f 2 \2n+1/_n+1 n+l* 


The load capacity ff for width b is given by 

kl 1 

W » 2b / p.dx + 2b / p-dx, 


(8 55 ) 
(8 54 ) 


which on using equations (S 55) and (8 54) gives 


where 


2bm^ 

^ ” rt^.O 


'^°o / 2n+1 „\n ^42 
n+2 ' 2n ^ ^ 


F - ( f )^"' (1 - 

2 e 2 


(8 55 ) 


(8 56 ) 


N= ^2 ^s 


The squeezing time t for reduo leg the film tnickness from an initial 

value h. to a final value h^ for a steady load W Is given by 
* dhg 

writing V = - as follows 
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_ 2 n +1 
“ 2 n 


2 bni 

[ 2— 

(n42)w 



^1 1 
i dh. 


(8 57) 


From equation (8 56} it con be very easily seen that both and 

8 

dF 

~ are negative for given h^i n and n ^is Implies that F 

decreases as b or k Increases 
8 


iDhue^ from equations (8 55 ) and (8 57 ) It can be concluded 

that both the load capacity Y/ and the time of approach t decrease 

with the inorease either In h^ In k !IhiB result is also obvious 

® h 

from the following approximate expression for t when ^ « 1 1 

f 


t 


2 ( 2 n+ 1 ) 
n 



1 



} 


- h k“^^ { T (8 58 ) 

2+1 1 + 1 

hj “ “ 

Ihe o1her obaraoteristios of W and t are similar to those obtained 
by ShUkla [1964b] 

B SQU^ZE FIIM OWO STEPPED PaRAII>SL CIPCULAR PLATES 

The squeezing flow of a power law fluid between two stepped 
parallel olrcular plates la shown In figure no (8 4 ) 

Hhe equation governing the flux In this case can be written 
(from equation (5 13) and putting b » 2 irr ) as follows 



155 


0 - 

2a+1 



j = 1 2 


( 8 59) 


where P-j ^2*^2 preseures and the fluxes In the 

regions h=:h^ 0 ;<x£1jH and h = hg, kR r ^ H respectively 

further, caa also be obtained from eo nation (5 15) as 


>r-^ = 2ir rV (e 60) 

ar 

AB in the previous c^se, since the flux is aero at r - 0 and 
Ql« at the step r s kS, we have on integrating equation (8 60) 

Qj = IT r^V (a 61) 

9!ben, from equations (6 59) and (8«6l) we get the equation 
deteimlning the pressure as fbllowsi 





1,2 


(8 62 ) 


Integrating this and using the boundazy conditions 


Pi e pg at r = KB 
Pg s 0 at r = R 

we have 


(8 63) 




HZ- (22111 

n+l 4n hj^ ^2 

(« 64 ) 
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n+1 


(2*^1 V)" (8 65) 

4n >9 


Obe load capacity 17 Is defined by 


Jrn D 

1 * / n, 2ir rdr + / p 2‘ifrdr , 

kP 2 


(8 66 ) 


udiich on using equations (8 64) and (8 65 )» gives 


^ ( igi . 


(8 67) 


where 


« - < * ( f ) (e 68) 


'2 s 


As In the case of parallel plates | the tLine of approach Is given 
by 


t . Isil. Z V <)hj 

4n ‘■(n+S)*''* ^ 2 


(8 69) 


Here the behaviour of & with respect to h„ or k is the same as 

that of F in the previous case Hence In this case also both the 

load capacity and the time of approach decrease as h or k increases 

!Qila result fbr tine of approach can also be concluded from the 

h 

8 

following approxlnate expreesion far t when "Jf” 
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t as 


r ^ p 

n ^Cn+3)»f 


” E: 


n r 1 
^n+1 1 + i. 

hf " 


1.1 

n 


} 


^-rhr 


2 + i 2 + i 
f 1 


->1 


(8 70) 


When h sa 0 la both the above OAsae the resulte are the seme as 

6 

obtained by Shukla [1964b] 

8 4 00NCU7SI0NS 

In this chapter the effects of Inertia and pressure on the 
oharaoterlstdne of squeeze film bearings wllh constant film thidmesB 
using power law lubricants are studied It Is found that the effect 
of Inertia Is to Increase the load capacity of the squeeze film 
bearing In both the oases On considering the linear varlatioh of 
consistency index with the pressure 1 e m - m^(l . ap)f It has 
further been pointed out that the load capacity Increases as the 
coefficient a of pressure Increases In the oase^ vdien Inertia Is 
neglectedi It Is seen that the squeezing tltoe Increases as a increases 

It has also been shown that In the absence of inertia^ and 
oonsistency variation with preesure, both the load capacity and time 
of approach decrease as the step height increases or as the step 
location moves away from the central position 



CHaPTJB - u 


OPTIMUM EXTERNALIil PRi!.SSimiSEO MA@fj:,TOH!a)ROI(XM<^lIC SEiffilMG 

As described la Chapter the problem of hydromagnetlcally 
lubricated externally pressurised bearing with constant film 
thickness has b en studied by Hughes and JIco [1i^62] f ^ukla and 
Prasad [1366 ] and others, Kamlyama [l969'-‘70-71 ] In this chapter, 

we consider the oase of an externally pressuilBed bearing vdth 
variable film thickness in the presence of a fon-unlfonn magnetic 
field applied in the axial direction The bearing surfaces are 
assumed to be non-conducting and no external current Is applied 
Due to external pres svudsat Ion, the lubricant flows symmstrioally 
in the radial dlreotlon and since the magnetic field is applied in 
the axial direction, a body force is generated in the film modifying 
the pressure of the lubricant The physical situation is illustrated 
in figure no (9 1 ) She calculus of variation technique, ae used 
by Rayleifh [1918] , Osteple and Young [ 1 962 ] and Valker and Osterle 
[1961 ] has been applied to determine the optimum profiles for film 
thickness and magnetic field for maximum load capacity of the 
bearing under a given supply pressure 

9 1 BASIC iQUATIOKSt 

Sie basic equation ^vernln^ the flow of a oonduotlng 
incompressible lubricant, under the usual as8i;)mptloa8 of hydro- 
magnetlo lubiicationi Shukla [196?] , Hugjhes and Jaloo [ 1962 ] , Is 




160 




3"U ir . 1 do 


u » — 


Bz' 


V dr 


(9 1) 


where Is functioa of r only and iii{r) s= h^B(r)(~)^ Is the 

Hartmann number 

Integrating equation (9 1 ) and usln«, the boundary conditions 


u = 0 at z - 0 

u = 0 at z = hCi^) 


(9 2) 


we get 


u » 


li- {cosh ^ - tonh I2L sinh ii - 
p<2 dr 2ho ^o 


1 } 


(9 3) 


abe flow fltix Q Is aeflned by 


r = / 2'irr udz ) 
0 


(9 4) 


which gives on usli^ equation (9 3) 




^3 ir ' 21^ 

Bearrangung equation (9 5) we get 


?h„ 


(9 5) 


dr 


yor 


1 


4irihJ ~ - tan*' ^ 
° Zi-V, 2ho 


(9 6) 
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It can be seen from the equation of continuity that is < 
with teepect to r Xbq bouudazy conditions for p are 

p ss p^ at r « 

P = Pa = ""e 

Integrating equation (9 6) and using the condition (9 7} we 


and 


Pa = fPi “ ’’a^ 


dx 


x(^ - tanh 


/ 


dx 


xcips. . tanh 


V 

1 a 


VO 




4irl 1 x(|!l - tann 

o 2 2 


dx - 

ry f 


where Y = Y(x)f x * ~ & b = — 
o ^1 ^1 


iiiquations (9 8) euid (9 9) can be written as 


and 




Pa ' (Pi - 


El 


/ 

I 


/ 

1 


FCi<,Y) 

X 

FQ V 

X 


dx 

— ] 
dx 


UQ 

W, 1 


iCLJ[Idx 


constant 


(9 7) 

get 


(9 8 ) 


(9 9) 


(9 10) 


- Pa 




(9 11) 
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Hie above expression gives the load oapaolty of the bearing for 
gl-ven fuactions Y(x) and mCx) 

Qbe following two cases are studied 

(a) Is presorlbed 

(b) Q Is presorlbed 

9 2 UAXXdOM lOaB CABACITf HW BYlCEBBAlOr ABH13D mSSORi^ IS 
IK&SCBIShD 

Here we wish to detexmloe the laagnetlc field and fl]Da thiokness 
pioflleB whlob maxiniae the load oapCK^ity of the bearing by using 
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calculus of -variations She following three cases are considered 

I when the film -tnlcKnese Is uniform and the applied magnetic 
field varies with radius, 

II when unlfoxm magnetic field Is applied and the film Ihlckaess 
varies with radius, 

III when both m-ignetlc field and film thickness vary with radius 

CASE I In this case, -the nondlmenslonal film -thickness Y is 
unity and hence from equa-fclon (9 14) the dimensionless load capacity- 
can be rewritten as 


b 

/ xFC! l)dx 

H i ^ , (9 15) 

>ra>,Pa)r? /” E£LiI dx > 

1 * 

where P, “ / x FC’ 1) dx dx (9 16) 

1 1 i 1 X 


and the function Iil(x) is to be determined for maximum W 
Now any change dM in M would make a corirespondlng change 6W 
in W Ibis change 6V can be written by using equa-tlon (9 13) 
follows 


where 


6W 


ff^(x) ts 


T > 
*1 1 


~ - Xf, 
X 1 


G,(x) — pa< 1) 6M dx 
* df 


(9 17) 


(9 18) 



1o4 


It oau te seen from ec^uatLon (9 12) that f(M|l) Is an Increasing 
function of M for n ^ 0 is and so P(j^, l) 21 0 11 ^ 0 

Qiis means that the sign of fV is given by the product of G-^Cx) 
and How the optiioim magnetic field is that value of 11 for 

which 5W la negative for any admissible 6 I 


Itollowln^ Rt-ylei^ [1918] we see that 6^(y) is positive 
for values of x close to imlty and negative fsr values close to ¥ 
Suppose the value of x at which &^(x) changes sign is then 

from equation (g 18) 



Ihe condition that dW is negative for any admissible 5M would 
be satisfied if we consider the following migcetlc field profllet 


d = 1 .1<x<a 
— — m 

ett <x<b 
o ' m — — 


(9 20) 


It can be seen that for 1 x dM oan only be positive and 

sinoe 0^(x) la positive In tills range, dW la negative Sbr the 
region a^ <, x <, b, dil can only be negative and since G^(x) la 
negative In this range dW Is again negative 3hus It has b en 
shovn that the optlmim profile fbr magnetLo field function i^ould 
be a step function 

Site load Capacity for this lype of mugnetio field profile la 
given by using equations {9 19) atd (9 19) as follows 
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W * 
m m 


(9 21 ) 


where a^ Is found by using equations (9 19 )« (9 20) as 


a. 


ft_ b 

/ xF(0 1) dx + / y F( * 1) dx 

2 'l 1 S 

In 


‘ /“" maid*,/ 

1 X 


pCq 1) 


dx 


n 


or 


PCIoD 2a;Jtna„ (aj l) 

F(0 1) “ . 2a2 ir 

m n a„ 
r' 


where P(M^i 1) and P(0|1) are the values of f( IjY) when Y ® 1> 

and M S' and M = 0 respectively Sttie equation (9 22) gives 

the relation between M and a for nuxlmum load cap^lty Since 

o za 

f(m ,1 ) increases as M increases, the function Sg also increases 
° a 

with The varl-tion of with is shown in figure no 

O do 

( 92 ) It is seen tiiat the function Sg inoreases -s increases 

for fixed b Since from equation (9 21 ), 13he optimvm load capacity 

2 

is proportional to a^ and therefore to for given b, it may be 

concluded that Wp, increases as increases 

Numerical example Ihe following example illustratas that the above 
chosen profile for 14 gives greater load capacity in caaparis on to 
uniform ma^mtlc field 


®(0,1) 


talcing the value of 10, in equation (9 22), we have 


10 41 = S„ 


Erom figure (9 2) we note that thrs vdlue 
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of Sg coairesponds tc 

ji = 3 16 -ind th€4i * 9 98 

m iQ IQ 

If we take ^ 0 throughout (or a constant m^^tic field is applied 
throughout ), then follovilng Hughes anu blco [19621 v/e have 


f/ = 


2Jln b 


5 


41 


Inus It Is concluded that V/ > V and the Increase is about 63% for 

IQ 

M = 10 
o 

CASH II In this oase uniform magnetic field is applied through- 
out the region 1 and Ihe film tnickness yCx) la a variable 

function of v te before, th*» expresoion for T can be \/rltten as 


b 

/ xFO^o Y) dx 

iLo.^ dl 

J 



(9 23 ) 


where 


and 



X P(!J^,y)dx 

P(.4 ,y) 

to 


(9 24 ) 


AS in Case I any change 6Y in Y will give rise to a corresponding 
change 6V in W 33160 from equation (9 23) given by 



b 


/ "oW 
1 “ 


tanh^c-^) 


6 


6Y dx 


(9 25 ) 
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$2 

where G2(x) = “ xf g 


(9 26) 


By applying, sliEdlar reasoolng as given in Case If It can be seen 
that the following profile for filjn xhicimess gL-yee the inaxiinum 
load capacity for uniform magnetic field 


Y = Y in 1<x<a 

o — — 

Y*1 in a<x<b 


(9 27) 


where Y is the ratio of the maximum to minimum film thloknass and 'a' 
o 

is the value of x at which GgCx) changes sign Ebr this profile also 
is negative for 1 ^ x ^ b and the maximum load is given by 


IL 


m 



£quatioas(9 23) , (9 27) and (9 28) give 


(9 2B) 



(9 29) 


and a is determined by 


PCMj 1) 2a^ tn a - (a^-1) 

"OoV'o.^ a2) - 2a2 »a i ’ 

a 


(9 30) 


ydiere the function ia Ihe seme as Sg where is repLaoed 
by a 
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Bie fUoctlon Sj is plotted with respect to and in figures 

( 92 ) and ( 94 ) Blgurea (^2) ^nd (9 4^ that the function 3^ increases 
as increases or appro^hes unity while figure (9 3) shov/s that this 

function decreases as M increases &r ^ 1 S^om this it is 

clear that t- increases as Y increases and decreases as 14 
bo 0 

increases Prom equation (9 28) it la seen that the optimum load 

2 a 

capacity is proportioral ta a and therefore to for fixed b, and 
so it may be conclvided that Increases as the step hel^t ratio Y^ 

increases and decreases as A Inore^ses for Y > 1 Thus > it may be 
remarked that a step beailng would give maximum load whmi there is no 
applied magnetic field in tl^ system 

CASE! Illi In the above two oases we have seen that for the maximum 
load capacity of tne bearing elthar the film thickness profile or the 
magnetic field profile should be a step finction Now the question 
naturally arises as to what happens when both non*-unifoxmity in magnetic 
field and step in film thickness are taken into account We investigate 
tbij question in the following* libr T->riable magnetic field m(x) and 
variable film thin mess y(x) the lo<^d capacity of the bearing can be 
rewritten as I see e4uatlon(9 14)] 


U 


/ xF(M,y)dx 


/ ££illdx 



(9 51 ) 


X 




ithere 


(9 52) 


b 

= j xI'CjAjY) ax 

1 

^ 1 

Now any change 6M In M and jY In Y will give net to a change 
aw in W Nron equation (9 51 )> 6V can be vnrltten as follows; 



(9 33) 
(9 34) 


Here also we batve to consider ihose vanations of jfl and Y such 
that the sign of dV la negative for maxlmuni W ^allowlz^ 

Haylelgh [1918] / the sign of Cr^Cx) la eeen to be poeiti-ve when x 
goes to 1 and It le negative when x approaches b 

let u8 consider the folio ing profllea for magnetio field m(x) and 
flln thiclcnasQ y(x) 

in = ^ , 1 < X < a 

— m 

M = jd^ a^<x<b (9 35) 

Y = Yq * 1 IX la 

Y = 1 , a < X < b 


where Is the wlue of X at which changes sign 

let us assume thot a It is seen that for the region 

1 1 X £ a» (^(x) ^0, ^ ^ negative which 
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implies that 6 A is ne^^itive Agaia, in the region 

since <^(x) ^ 0| ~ >, 0 for il L Of < Of 6iS. < 0 and 5y is 

poBitiv6f so that chc sign of ^ i is clearly negative Xiet us now 

examine the sltaO of si/ in the region a <,x ^ In this range 

/ \ ftp dP 

GjVx) is positive, “ is positive, &U is positive, is 

negative and 6Y is positive, so the sign of 6W may be positive 

or negative [sec e<^uatioa (9 33)1 Since for maxi>suBi W, 6 W is 

to be negative in tiiis region also, we should take a situation which 

makes > 0 and 5Y < 0 Qlhia is clearly satisfied if we make 

a tend to a Ihe same argument holds in case < a ^Qius the 

above chosen profiles give 6W < 0 tnrou^out the region, then 

a B a 
m 

Since G^(x) = 0 at x « a^ » a, from equations (9 31) and 
(9 34) we have I 



®3 2 
^ = a 

3 


(9 36 ) 


How, from equations (b 32), (9 35) and (i# 36) we have the expression 
determining a as follows 


j, _ , (a^-l) p(0,Yq) •<- J) P(M^,1) 

^ P(0,Y^)Ana+ p(Mp,l) tn t 


or, 


2a"ltna..(a^.l) 


(9 37) 


(9 38) 
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It can be again noted here, that the function — Increases as 

Ho,y^) 

or Y or both increase Funher, as S. increases as 'a 'increases 
for given b, it is concluded | as befc»?e that increases as 
or or botn increase It is remarked that this maximuai load is greater 
than that of the previous two cases | and thus it is desirable to have 
a magnetohydrodynamic externally pressurised bearing with stepped film 
thicKneas and the magnetic field snould be unifomOy applied only In 
the minimum film thickness region [see figure no (9 5)] 

9 3 lOAD xHM THh MOW RASE. Q 18 IR±.SCBlBJDt 


In tnls case the load ca|.aclty V is written from equations 
(9 11) and (9 U) as 


yO b X ^(x) 

“■"Jiir { ' w xa . t„n , t'WVCx) 

0 2 2 


(9t59) 


Since F(11|Y) increases as ft increases for given it m£Qr be 
infeirred from equations (9 11) €jnd (9 39) that for a given both 
p. - p and W increase with increase of M It is further observed 

1 E 

njEV 

that Bioce - tahh — lacreases as y IncreaseB for fixed lly W 
increases as Y decreases i e the step height ratio approached unily 
ThuBi for a given flow flUKy step in the film thickness is not desirable 
(Oils result is true even when there is no magnetic field Sbr Y = 1 
(constant) we obtain all the results from here as investigated by 
Hui^ee and Blco [1962] 
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9 4 CONCIIJoIONS 

In tbe ll^t of the above analysis th^ following conclusions 
may be drawn 

(l) At a constant applied pressure, Ihe magnetic field should be 
nonunl orm for optimum load c<^pnclty if the film thickness Is constant 
This maxlnum load Increases as the strength of the negnetic field 
Increases 

(•^) Bor uniformly applied stagnetic field the film thickness 

should be a step function If the load capacity Is to be maximum and 

this maximum loud capacity Increases as the step height ratio 

Increases and it decreases as M Increases for Y > 1 Thus, for 

0 o' 

maxlifium load capacity it is desirable not to apply any magnetic field 
in the case of step bearing 

(?) It Is seen In Cass IH that, for tne inaxlmum load capacity of 

the bearing, both the film thickness and the m.gn«:tic field functions 

should be step functions huvlng th same step location [ see fig 

ao (9 5 )] Ibis maximum load increases as H or Y_ or both 
' o o 

Increase 

( 4 ) ^en the flow rate of the lubricant Is constant, both the 
recess pressure and load capacity Increases with increase In the 
strength of the magnetic field Also the load capacity decreases 


as the film tnlcknees ratio Increases 



CHAPTIH ~ X 


MAGlfE!K)Hn«Oiflf^JUlIC EXPERNAUaf PR2SSUBISED POHOUS THEHTST BBARUTO 


In Chapter IX t the magpetohydrodynanilo externally pressurized 
bearlQg has been studied when extexn:il pressure is applied at 3 
Central recess In this chapt r^we study the cnaructexistlcs of 
magnetofaydrodynnmlo thrust bearings vdiere one of the surfaces is 
porous Por the s^lce of 8ii%llclty the porous matrix has been 
assumed to consist of a system of capillaries of very small radix 
which restrict the lubricant to flow only In the direction normal 
to the plates !Qus assumption seems to be reasoruable due to recent 
advances In cernmlc technology, Hslng [1971] A constant magnetic 
field has betn applied throu^out the bearing system In the transverse 
direction vdilch modifies the pressure In the film Since the applied 
magnetic field le in the some direction os the flow of the lubricant 
in the porous matrix, the velocity of the lubricant in the matrix 
remains unaffected and Is still to be governed by the usual Darcy's 
law ^Hslng [1971] 

Here we study the following two cases of MHD pomus thrust 
bearings 

(l) lUHD porous thrust Infinite plates, 

(lx) lifilD poroios bhrust circular plates 

Sitollar noomagnetlo Ctoses have been studied by dslng [ 1971 ] } aad 
Mori et al [1965] 
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10 1 15HD lOROUS THHUSI ILaIxiSs 


Consider the How of an Incompressible lubricant between two 
Insulating iK.ra.ll 1 plates In the presence of a transrers ly applied 
magnetic field The upper plate Is porous through which the lubrlc<mt 
Is pressurised Into the bearing and thus generating a symmexrlc 1 
pressure dls rlbutlon in the film Ihe physical situ.^tion is shown 
In figure (lO l) 


Ihe basic equ..tlon governing the flow of the conducting lubricant 
under the usual hydromjgnetlc assumptions ^ Kunna [1965] , Is 

2 


a^u 

ay2 

1 

ra,T 


h2 


u » , 

y <?x 


where M — £b C^)** Is the Hartmann number 
for the velocity are 


(10 1 ) 

{{he boundary conditions 


u « o ^t y » 0 
u = 0 at y = h 


(10 2 ) 


Solving equation (l0 l) and using (lo 2) we get 


u 



[{rosf 


- U - £2211 
h in 




(10 3) 


The flov flux, «i for width b is defined by 

h 

a b / udy 
0 


I 


(10 4) 






kB 1 
iPe 

IpJ 
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FIG 10 1 MHD- POROUS THRUST PLATES 
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which on using equution (lO 3 } gives 

3 
^r: 

V 


dy 2 


(10 5) 


Ghe equation of continuity is given by 


0 

3x 3y 


(10 6 ) 


Integrating equation (10 6) and using (10 4) and the conditions f 


V « 0 

at 

y = 0 

Y a -V* 

at 

y = h 


(10 7) 


where r' is the velocity in the porous mesh) we get 


h 


(10 8 ) 


From equations(lO 5 ) and (IO 8) we get 

(’“S) 

She yeloollgr v' in the porous mesh is given by Bazoy’s law as 
[Hsing (1971)] , 


at _ IJ! 


p H 


(10,10) 


vhere ^ is the permeability) u is the viscosity of the lubricant 
and H is the thleloiess of the porous matrix equations (IO 9) 

and (10 10) we obtain 
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Solving equnxlon (10 11 ) and using the conditions 


we get 


1^-0 at X = J 
dx 


p = 0 at X = -^ I 


cosh 


2Bx 


V - D^Cl 


CC •" f 


-) , 


vdiez^ 


8?. « 


^ - tjoh j 


fflib load rapacity / of tho bearing for width b ie given by 


V2 

W = 2b / pdx 
U 


which on using equation (l0 13) gives 


\f__ ^ . tanh g 

blp^ ' “ " 


Since 


d tanh 6 

de ^ 0 


sinh 20 > 20 

**0^ cow.h^0 


(10 11 ) 


(10 12 ) 


(10 13) 


(10 14) 


(10 15) 


(10 16 ) 
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tsuih 3 

Is positive for a > 0 1 therefore 1 itioreases as 3 

p 

increases Further from equation (l0 14) it can be seen^ as 

increases as M increases, that 3 increases as M 


f -tanhf 


or 4 increases Xhus it can be concluded from equation (l0 15) that 
W increases -.8 ^ or (j> or both increase 


By putting M = 0 in equation (l0 15) we gut the following expressions 
for the load Ofe.p-.city for the non-^agnetlc case which is the same as 
obtained by Hsing [1971] when inertia effects are neglected in his 
Ct^se 

w = 1 - , (10 17) 

6 = =■ ^ 

° Hh^ 

For 11 = 10 and 3^= 5 ) the load increase from the non-magnetic oaee 
is about 


1 0 2 SQUEEZING JfePFECT 


Now consMer the above case when the upper plate is having 
a relative nozml velocity V due to loading, then the conditions 
(10 7) are modified to 


V = 0 at y « o 

V ss -V-v' at y = h 


(10 18) 


!l!hun following the same procedure as before the expressions for 
pressure and load capacity can he written as, 
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" = a - 


cos 3 


(10 19) 


- . V , Cl *^1(1 


(10 20 ) 


where = •— , V and (p^+ ^ ) (l- ) bL 

Now from equai:loa (10 20 ) It oen be seep th£>.t loud capacity locreaeeB 
as V or a Increases Qbe squeezing time t for reducing the 
film thickness fiom to la given by 


1 


(10 21 ) 


3>t8nh 3 


“1 — h 

where h. « xr~ h = t— 

1 hg hg 

Since 1 - is an inoreaalng function of 3 it can, be seoi 

from equation (l0 21 ) that t Inoreasee as M increases for a given 
When ^ = 0} the expressions for W and t are the seme 
as obtained by Sbukla [196^] for large and small Hartmann numbers 

30 see the effect of ^ only) we consider the case when M » 0 
In this case the expressions for W for small ^ is given by 


V I fv + (l --^K) 5] 


(10 22 ) 


where 



from equation (10 21 ), the approximate 'value of t in the form of an 
Integral is given by 


t 


u 






(if' 23) 


equations (l0 22) and (10 23) it may be concluded that W 

increases as ^ increases provided V < 3/6 m - which may be 

VL 

satlafied by a suitable choice of 

Jigalntfrom equation (l0 23) it may be inferred that t increases 


as ^ increases prcvided W^< -g- When M is very large ^ from eqviations 


(lO H) and (10 2l) j«b hsT«3, t * which implies th-t t decreases 

^ Increases 

10 3 MHD POROUS IHRUST OIRCULaB PLaSCtS 

In tnls case also the upper plate is porous jiad the lujrlcant 
is pressiucised into the beanng which generates a symmetrical pressure 
distribucion She physical configuration is illustrated in ftgure no 


(10 2 ) 


Ohe bas3c equation governing 'the flow of the incompressible 
lubricant wider 'the lubrication assumptions, Hu^es and Elco [1962] , 
is given by 


8^u 1 ^ 


(10 24) 


where M is the Eortmann number, u is the -velooity in the radial 
direction and ^ is a function of r only 3he boundary opodltions 
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for th« velocity u ar 

u = 0 at z - 0 

(10 25) 

u = 0 at a = h 

Solvliig equation (10 24-) with the help of the conditions ( 10 25) 
we obtain 


h* do , .Hu , 
VM 


The flow flux Q le defined by 


Cosh M-1 
slnh id 


sinh 



h 

0 = / 2irr udz 
0 


which on using equation (10 26) gives 
3 

^ 4Trrh dp M Mn 

Q ^ [tanh ~ - ~3 

The equation of continuity In this ease Is given by 


(10 26) 


(10 27) 


(10 28) 


i|^Cou,.|i=c (10 29) 

Integrating equation (10 29) with respuct to a and using equation 
(10 27) and the conditions 

V * 0 at a = 0 

(10 30) 

V = -V* at 3 = h 


we get 


^ « ?irrv» 

ax 


(10 51) 
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wh(.re v' is the velocity of tht. lubricant in the porous matrix 
from equations (l0 28) and (l0 31 ) we get 

Air? ^ f i1 Ht> 

^ Ctanh ~ -) ^ (r . '’vr v* , (lO 52) 

yf j 2 A ar dr 

13ie velQoity y* in the porous matrix Is again ^ given by 

^ P “ P 

(10 33) 

from equaTions (IO 32) and (IO 3:>) we get tae final equation 
dete rminin g the pressure as 



f R d f dp. 2vr ^Pe 
2 “ 2^^ dr dr^ “ y V 


(10 34) 


Equation (10 34) can be written as 


.2 d-^P . ^ dP 


«?? 




(10 35) 


where E » p^~ p, ^ ^ 


2Hh-' f - tcnh I 


and f » or 


Ihe boundary conditions for p are 


=* 0 at r = 0 

p ss 0 at r = R 


(10 36) 


OJ7 equlval ntly, 


^ » 0 

dr 


at 6 « 0 


I = p^ at 6 - OR 

t3 


(10 36)' 
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The solution of equation (10 3^) can be vFilttsn as 


P = A I^(C) + (1037) 

vAiere A and B are constants and are modified 

Bessel function of ord r zero of first and second kind respectively 
Since the pressure la finite v.t 5 = 0 ^nd ^ infinite 

singularity at ?= 0, we get B = 0 Then from equation (l0 37) we have 

p = aIq(c) (1038) 


iifter determining the constant A from the second boundory oondltion 
(10 36) the expression fbr pressure p can be written as 


V = 


loCO 

ilw ^ * 


(10 39) 


where s on 


She load capacity W is given by 


VJ a / 2 t(t pdr 

r 


(10 40) 


From equations (IO 39) and (10 40) ve get 


1/ = 


= [1 




irn^P" 




3 


She variation of the function 1 - 


2 ^^ ( 6 ^) 


(10 41) 


is ehowi in figure no 


(10 ^)ffTom vhlch it can be seon that this function increases as 
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2 6R 

3. incre>cs«s Since 6- « — sr — ^ Increase as A or 

' 1 2St? f-tenh| 

M or both Increase, -cherefore it can be concluded from equation 

(10 41 ) that W Increases as <|) or M or both increase 


10 4 SQU^ZllfG *a!’?iaOa? 

We now consider the above case when there is a relative normal 
Vdloolty V between the plates Ihe conditions (l0 30) are now 
modified to 


V » 0 at z = 0 

(10 42) 

V -v-v' at z = h 

frooeedlng as before the e:(pre88ion8 for pressure and load capacity 
are given by 


o 




WHV, ,, 

+ -7-) (1 


loCC) 
-= — ] 
IpCSj) 


and 


V 




2 Ii(3i3 




-] 


(10 43) 


(10 44) 


making dimenaLonless, we hove from equation (10 44) 

'v V 2 Ii(3i) 

« ^■= (1 + r ) [1 •] ( 1 «> 45 ) 

It is seen fxom equation (l0 45) as before that the load capacity 
increases as V or M increases 
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!Qiii squeezing time t for reduoing the film thickness from h^ to 
hg is given by 


« I A. 

p’lo ^ 1 


(10 46) 


?llCPl) 


where funotions of h Sixice 

1 - ^ is an inoreaslng function of (see fig 10 3) it la 

seen from equation (10 46) th£.t t increases as M Increases for a 
given 4 By taking 4 =0 in the above expressions for W and t 
we get all the results as obtained by Shukla [1965] for large and small 
Ifortmann numbers 

In order to see the effect of 4, the corresponding load 
W_ for small and M can be written as follows 

V 1 


~ [V + (1 ) 4] 


(10 47) 


where = K,4 and K, =7 

^ ’ 2h5 f-tanhf 


It ma^ be seen from 


equation (10 47) that Inoreasss as 4 increases provided 


12h^p — tanh ^ 

*^e 2 _ 2 


suitably adjusting the applied pressuro p 


and this oondltloh. may be satisfied by 
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Again for email t Is given from equation (10 46) as 



1 


dl 




(10 48 ) 


Since > 1 ^the IntegranSl and hence the integral increases as ^ 

increases provided ^ ^^<3 Inqjlies that the time of approach 

~ - 3 

is delayed due to increase in 4> provided W , < y 

V 4 

Further for large M and hence from equation (10 46), t 
can be approximated to 

t =1 ( £i- 1) 

* *v- 1 

and decreases as ^ Increases fcr a given 


10 3 GONCBJSlONSs 

Ih this chapter, we have studied the effects of transversely 
applied magnetlo field on porous thrust bearings by ooneiderlng both 
rectangular and circular configuraticns It has b an shown in these 
cases that the load capacity increases as supply pressure, pemeability 
and the strength of the applied magnetic field increase Xbe effects 
of squeezing have also been studied in the above mentioned oases and 
it has been pointed out that the load oapaoity Inorsases as squeeze 
veloolty Increases and this Inorease is enhanced as the strength of 
the applied magnetlo field increases fbr a fixed ^ It has further 
been pointed out that the time of e^proaob also Increases as *016 
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strength of the applied magnotlo field mcreases for a prescribed 
In the case when squeezing Is taking placc^ the load capacity 
Increases as (j) increases provided certain restrictions are placed 
on V as htjB b«en pointed out earlier Slnn-larly the time nf 
approach also Increases as <(• increases for a given M subject to 
certain oondltions on 
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